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3.1  The  criterion  of  stability. 

Our  discussion  of  the  general  theory  of  stability  of 
equilibrium  will  be  based  on  the  energy  criterion  for  conserva¬ 
tive  systems.  As  we  have  already  seen  in  section  1.3^  the 
implied  existence  of  a  potential  energy  is  not  too  severe  a 
restriction  in  the  case  of  elastic  solids  or  structures.  An 
internal  potential  energy  always  exists  in  isothermal  or  adiabatic 
deformations  of  an  elastic  solid,  and  the  external  loads  on  a 
structure  are  also  conservative  for  most  static  loading  condi¬ 
tions  encountered  in  practice. 

In  the  present  chapter  (v/hlch  is  largely  based  on  the 
analysis  developed  in  [10,  ch.2])  we  shall  confine  our  attention 
to  an  investigation  of  the  stability  of  a  particular  configuration 
of  equilibrium  under  the  action  of  completely  specified  external 
loads .  V/e  shall  denote  this  supposedly  knov/n  equilibrium  con¬ 
figuration  as  the  fundamental  state  or  state  I .  The  energy 
criterion  of  stability  of  this  fundamental  state  may  then  be 
expressed  in  somewhat  vague  general  terms  by  the  statement  that 
the  existence  of  a  proper  relative  minimum  of  the  potential 
energy  in  state  I  is  a  necessary  and  sufficient  condition  for 
the  stability  of  this  configuration.  We  shall  nov;  examine  a  more 
precise  statement  of  this  energy  criterion  which  is  appropriate 
for  its  application  to  continuous  systems. 

The  investigation  of  the  minimum  properties  of  the 
potential  energy  in  state  I  requires  its  comparison  with  the 
potential  energy  in  all  possible  adjacent  states  II  within  a 
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certain  neighbourhood  of  the  fundamental  state ^  to  be  defined 
presently.  The  possible  adjacent  states  are  characterized  by 
the  property  that  they  satisfy  the  geometric  conditions  of  support 
for  our  system.  The  displacement  field  u  which  carries  our 
system  from  state  I  to  a  possible  state  II  v/ill  be  called 
kinematically  admissible.  The  increment  in  potential  energy  of 
the  system  in  the  transition  from  state  I  to  state  II 

PlI-Pl  =  P[u]  (3.1.1) 

is  a  functional  of  the  kinematically  admissible  displacement 
field  u.  It  consists  of  a  line  integral^  a  surface  Integral,  a 
volume  integral  or  a  sum  of  such  integrals,  whose  integrands  all 
depend  on  the  local  displacement  vector  and  its  spatial  deriva¬ 
tives.  We  shall  say  that  a  displacement  component  or  derivative 
appears  explicitly  in  (3.1.l);»  if  at  least  one  of  the  integrands 
depends  explicitly  on  this  component  or  derivative.  Let 
u^(i=l,2,3)  denote  the  continuous  Cartesian  displacement  compo¬ 
nents,  and  u.  ,,  u.  . (i, J,k=l,2, 3)  their  continuous  partial 

derivatives  v/ith  respect  to  Cartesian  coordinates .  We  define  a 
neighbourhood  f  .  ..  ]  of  the  fundamental  state  I  by  the 

set  of  possible  states  II  subject  to  the  condition  that  the  compo¬ 
nents  and  derivatives  of  the  kinematically  admissible  displacement 
field  from  state  I  to  state  II  which  appear  explicitly  in  (3.1.1) 
satisfy  the  inequalities 

I ±  e  ,  I is'  ■>  I j J  £  s"  ’  •  ••’ 


(3.1.2) 
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where  g,  gS  &" j  • • •  are  positive  constants.  A  necessary  and 
sufficient  condition  for  the  stability  of  equilibrium  in  our 
fundamental  state  I  is  nov/  the  existence  of  a  non -vanishing 
neighbourhood  f gj  g*,  of  state  I  v/ith  the  property  that 

the  energy  increase  (3.I.I)  is  non-negative  for  all  possible 
states  II  in  this  neighbourhood.  In  other  words ^ 

P[u]  2  0  (3.1.3) 

for  all  kinematically  admissible  displacement  fields  which 
satisfy  (3.1.2)  for  some  set  of  positive  constants  g,  g’^  g" , • • • , 
is  a  sufficient  condition  for  stability.  On  the  other  hand,  if 
in  any  neighbourhood  {gjgSg”j...}  ^  no  matter  how  small  we 
choose  the  constants  gjgSg":>...j  an  admissible  displacement 
field  u  exists  for  which  P[u]  is  negative,  then  the  equilibrium 
in  the  fundamental  state  is  unstable . 

The  reader  who  is  familiar  with  the  calculus  of  varia¬ 
tions  will  note  that  our  formulation  of  the  minimum  property  of 
the  potential  energy  in  a  stable  configuration  of  equilibrium  is 
equivalent  to  the  definition  of  a  so-called  weak  minimum  in  the 
calculus  of  variations  [e.g.  3,9]-  We  have  two  reasons  for  our 
preference  for  the  definition  of  a  weak  minimum  over  the  more 
exacting  requirement  of  a  strong  minimum.  The  latter  definition 
would  restrict  a  neighbourhood  of  state  I  only  by  imposing 
bounds  on  the  displacement  components  u^  themselves,  and  not  on 
their  derivatives .  Prom  a  physical  point  of  view,  however,  the 
boundedness  of  the  derivatives,  in  so  far  as  they  appear 
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explicitly  in  the  energy  increase  (3-1. l)j  is  equally  necessary. 
This  physical  argument  is  seconded  by  the  fact  that  the  mathe¬ 
matical  theory  of  strong  extrema^  in  particular  in  the  case  of 
multiple  integrals,  is  far  less  well-developed  than  the  theory 
of  weak  extrema.  A  detailed  discussion  of  stability  conditions 
on  the  basis  of  a  requirement  of  a  strong  minimum  of  the  potential 
energy  would  present  nearly  insurmountable  difficulties  in  the 
present  state  of  knowledge  in  the  calculus  of  variations.  It  is 
indeed  fortunate  that  the  physical  requirements  of  stability 
permit  us  to  restrict  our  investigation  to  the  conditions  for  a 
weak  minimum. 

The  particular  form  of  expression  (3.1.1)  for  the 
energy  increase  in  the  transition  from  the  fundamental  state  I 
to  some  possible  state  II  is  immaterial  for  our  purposes  in  the 
present  chapter.  In  fact,  we  shall  deal  here  with  some  basic 
problems  in  the  calculus  of  variations  in  a  more  general  context, 
although  we  shall  employ  the  terminology  of  stability  theory. 
Nevertheless,  a  better  understanding  of  the  theory  may  be 
facilitated,  if  we  mention  at  least  some  explicit  examples  of 
energy  expressions  in  order  to  illustrate  the  otherwise  somewhat 
abstract  basic  theory.  In  this  connection  it  will  also  be 
convenient  to  anticipate  in  appropriate  instances  some  results 
of  a  later  detailed  discussion  of  these  examples. 

As  a  first  example  we  take  the  well-known  buckling 
problem  of  a  simply-supported  straight  bar  of  constant  cross- 
section  under  the  action  of  central  compressive  end  loads  N 
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(fig.  3.1).  Tlie  fundamental  state  I 


Pig.  3.1 

Buckling  of  bar  in  compression 

which  is  the  subject  of  our  investigation  of  stability^  is  the 
straight,  undeflected  configuration.  We  assume  the  center  line 
to  be  inextensible .  A  material  point  on  this  center  line  is 
identified  by  its  distance  x  to  one  end  in  state  I,  which 
distance  equals  the  arc  length  in  the  deflected  state  II.  We 
confine  our  attention  to  deflections  w(x)  in  one  of  the  principal 
planes  of  the  bar.  The  elastic  energy  is  zero  in  state  I.  It 

o 

amounts  to  ^  Bv.  per  unit  length  of  the  center  line  in  state  II, 
where  h  Is  the  curvature  and  B  the  flexural  rigidity.  The  energy 
of  the  prescribed  end  loads  equals  the  product  of  N  and  the 
distance  between  the  ends  of  the  bar.  Hence  we  obtain  in  this 
simple  case  expression  (3.I.I)  for  the  energy  increase  in  the 
form  of  a  single  integral 
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=  P[w]  =  B  +  N(/iI7'2  -  i)]dx  ,  (3.1.4) 

where  primes  denote  differentiations  with  respect  to  x.  In 
order  that  our  deflection  field  be  kinematically  admissible,  we 
must  have  w(0)  =  w(i)  =  0.  Equilibrium  in  the  fundamental  state 
I  is  stable,  if  and  only  if  positive  constants  g*  and  g"  exist 
such  that  for  all  kinematically  admissible  deflection  fields 
w(x)  v/hlch  satisfy  the  inequalities 

Iw'l  I  g'  .  Iw"l  <  g"  ,  (3.1.5) 

the  energy  increase  (3.1.^)  is  non -negative 

P[w]  >  0  .  (3.1.6) 

As  our  second  example  we  choose  the  three-dimensional 
stability  problem  of  an  elastic  body  in  a  state  I  of  initial 
stress  described  by  a  symmetric  tensor  (i,j=l,2,3)  referred 
to  Cartesian  coordinates  x^  (i=l,2,3) .  We  denote  the  vector  of 
body  forces  per  unit  volume  in  state  I  by  X,  its  components  by 
X^,  and  the  vector  of  surface  tractions  on  the  part  of  the  surface 
Sp  where  the  loads  are  specified  by  £,  its  components  by  p^ .  We 
shall  confine  our  attention  to  the  case  of  dead  loading,  in  which 
the  prescribed  loads  on  a  material  element  of  the  body,  specified 
loy  the  vectors  X  dv  and  £  dS,  do  not  change  in  either  magnitude 
or  direction  in  the  transition  of  the  body  from  state  I  to  some 
other  possible  configuration  II.  If  the  surface  displacements 
are  prescribed  on  a  part  of  the  surface,  a  kinematically 
admissible  displacement  field  u  with  components  u^^  must  vanish 
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on  S^.  The  (additional)  deformation  of  the  body  in  any  kinemat¬ 
ically  admissible  displacement  field  from  state  I  to  state  II  is 
described  by  the  strain  tensor 


(3.1.7) 


where  the  siommatlon  convention  has  been  employed  for  repeated 
subscripts .  An  appropriate  expression  for  the  increase  in  energy 
(3.1.1)  will  be  derived  in  chapter  6  for  small  (although  not 
infinitesimal)  deformations  of  a  solid  which  obeys  Hooke » s  law . 

Tlie  resulting  expression  reads 

p^u^dS  , 

P 

(3.1.8) 

where  Q  is  the  shear  modulus  and  v  is  Poisson *s  ratio.  Equilib¬ 
rium  in  the  fundamental  state  I  is  stable ^  if  and  only  if  positive 
constants  g  and  g*  exist  such  that  the  energy  increase  (3.1.8) 
is  non-negative  for  all  kinematically  admissible  displacement 
fields  u  which  satisfy  the  inequalities 

luj^l  Is  ,  I  Is'  (3.1.9) 

everywhere  in  the  body  and  on  its  surface. 


■[u]  =  J 


Ky 


+ 


l-2v'^^hh 


(Yv 


^)2}  -  X^u^]dv  -  J 
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з. 2  Some  necessary  conditions  for  stability. 

In  order  to  apply  our  general  energy  criterion  we 
assume  that  the  integrands  of  all  Integrals  appearing  in  the 
energy  Increase  (3.1.1)  possess  in  a  neighbourhood  (3.1.2) 
continuous  partial  derivatives  with  respect  to  all  arguments 

и .  ,  u  .;/  u.  etc.,  which  appear  explicitly.  For  the  sake  of 

1  1^  J  ^3 

brevity  we  shall  denote  these  arguments  by  y^(X=l,2,  . . . ,q) ,  and 
their  aggregate  by  where  q  is  the  number  of  argimients  u^,  u^  ^ 
etc.  which  appear  explicitly  in  (3.1.1).  Let  F(^)  be  a  typical 
integrand.  We  denote  partial  differentiation  v;ith  respect  to  an 
argument  y^  by  a  subscript  X,  preceded  by  a  comma,  e.g. 


P 


X 


etc . 


(3.2.1) 


If  we  employ  again  the  summation  convention,  now  from  1  to  q  for 
a  repeated  Greek  index,  we  may  write  the  Taylor -exp ansi on  of  the 
typical  Integrand  in  the  form 


f(£)  =  y 


.  1__  ^1 
m!  ^ 


.yS 


(3.2.2) 


where  the  arguments  q  of  the  m-th  order  derivatives  are  some 
positive  fractions  of  the  argimients  y^.  We  collect  the  integrals 
whose  Integrands  are  homogeneous  functions  of  their  arguments  of 
the  same  degree  n,  and  we  call  the  corresponding  sum  of  integrals 
this  way  we  obtain  the  so-called  Taylor-expanslon  of 
our  energy  Increase  functional  (3.1.1)  in  the  form 
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P[u]  “  ^2^—^  +  ...  +  ^  (3.2.3) 

where  the  asterisk  attached  to  the  last  term  Indicates  that  it 
is  the  m-th  order  remainder. 

It  will  be  convenient  to  assume  also  that  all  geometric 
conditions  which  restrict  the  class  of  kinematically  admissible 
displacement  fields  u,  are  linear  and  homogeneous  in  the  displace¬ 
ment  components  u^  and  their  derivatives  .  Any 

linear  combination  of  kinematically  admissible  displacement 
fields  ia  then  again  such  a  field.  In  particular,  we  consider 
the  linear  one-parameter  family  of  displacements  '  au,  where  u 
is  some  kinematically  admissible  displacement  field,  and  a  is  a 
parameter  independent  of  the  coordinates.  A  necessary  condition 
for  stability  is  now  the  existence,  for  every  kinematically 
admissible  displacement  field  u,  of  an  associated  positive  number 
k  such  that  the  inequality  !a|  k  ensures  the  inequality 

+  PgCau]  +  ...  + 

-  csPj^Iu]  +  a^Pgtu]  +  ...  +  0(a"*)  g  0  .  (3.2.4) 

IT) - - 

'A  separate  investigation  is  required  in  cases  in  vdiich  the 
geometric  conditions  are  nonlinear.  We  may  then  consider  a 
nonlinear  family  of  admissible  displacement  fields  u(a)  which 
reduces  to  zero  fora=0,  but  we  sheill  not  pursue  the  investiga¬ 
tion  of  these  exceptional  cases. 
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We  employ  the  symbol  O(a^)  here  and  in  the  sequel  for  any  quantity 
which  tends  to  zero  as  for  a  — >  0.  In  the  present  instance  it 
expresses  the  fact 

|a“^P*[au]|  <00  for  a  0  ,  (3.2.5) 

which  is  an  immediate  consequence  of  the  boundedness  of  the  m-th 
order  partial  derivatives  of  the  typical  integrand  P(^)  in  the 
neighbourhood  (3.1.2). 

The  first  necessary  condition  which  follows  from  (3.2.4) 
is 

P^Eu]  =  0  (3.2.6) 

for  every  kinematically  admissible  displacement  field  u.  This 

condition  expresses  the  principle  of  virtual  work  for  the  equilib- 
ri\mi  of  our  conservative  system  in  its  fundamental  state.  Since 
we  have  already  assumed  that  state  I  is  a  configuration  of  equi- 
librium,  eq.(3.2.6)  does  not  provide  any  new  information. 

The  second  necessary  condition  for  stability  obtained 
from  (3.2.4)  is 

2  0  •  (3.2.7) 

It  may  be  expressed  in  the  form  that  the  second  variation  of  the 

potential  energy  in  state  I  must  be  non-negative  for  every  kine¬ 
matically  admissible  displacement  field. 

If  the  necessary  condition  for  stability  (3.2.7)  is 
satisfied,  we  have  to  distinguish  between  two  cases .  In  the 
first  case  the  equality  sign  in  (3.2.7)  holds  only  for  an 
Identically  vanishing  displacement  field  u.  The  second  variation 
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of  the  energy  is  then  positive  for  every  non-vanishing  admissible 
displacement  field,  and  this  second  variation  will  be  called 
positive  definite.  In  this  case  we  cannot  infer  from  (3.2.4) 
any  further  necessary  conditions  for  stability.  The  second  case 
arises,  if  the  second  variation  is  again  always  non-negative  but 
if  it  takes  the  value  zero  for  some  non-vanishing  displacement 
field,  say  u^.  The  second  variation  is  called  positive  semi- 
deflnlte  in  this  case,  and  we  obtain  two  further  necessary  condi¬ 
tions  for  stability  from  (3.2.4)  in  the  form 

=  0  ,  (3.2.8) 

>  0  .  (3.2.9) 

Equation  (3.2.8)  and  inequality  (3.2.9)  must  hold  for  every 
kinematically  admissible  displacement  field  for  which  the 
second  variation  vanishes . 

If  the  additional  conditions  (3.2.8)  and  (3.2.9)  in 
the  case  of  a  semi -definite  second  variation  are  both  satisfied, 
we  have  again  to  distinguish  between  the  case  in  which  the 
fourth  variation  P|^[u^]  is  positive  definite  (i.e.  positive  for 
all  non-vanishing  displacement  fields  u^  for  which  PgLu^]  =  O), 
and  the  case  in  which  the  fourth  variation  is  semi-definite .  In 
the  first  case  we  obtain  no  new  necessary  conditions  from  (3.2.4), 
in  the  second  case  we  infer  additional  necessary  conditions  for 
the  fifth  and  sixth  variations,  similar  to  (3.2.8)  and  (3.2.9)^ 


etc . 
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V/e  emphasize  again  that  our  discussion  up  to  now  has 
resulted  only  in  a  nuraber  of  necessary  conditions  for  stability. 
The  question  to  what  extent  these  conditions  are  perhaps  also 
sufficient  conditions  for  stability  is  still  completely  open  at 
this  stage;  this  question  will  be  discussed  in  detail  in  the  next 
sections .  In  the  conventional  treatment  of  stability  theory  on 
the  basis  of  the  energy  criterion  the  discussion  is  usually  not 
pursued  beyond  (3.2.7).  It  is  then  taken  for  granted  that  a 
positive  definite  second  variation  of  the  energy  is  also  a 
sufficient  condition  for  stability.  This  complacent  view  is 
presimiably  based  on  a  tacit  assumption  that  (3.2.4)  also  supplies 
sufficient  conditions  for  stability.  This  tacit  assumption, 
however,  is  gravely  in  error.  We  shall  see  later,  in  section  3.6, 
that  equilibrium  in  the  fundamental  state  may  very  well  be 
unstable,  in  spite  of  the  satisfaction  of  (3.2.4)  for  every 
kinematically  admissible  displacement  field.  For  example,  even 
in  the  case  of  a  positive  definite  fourth  variation  P]|[u^]^ 

(3.2.8)  and  (3.2,9)  are  not  sufficient  conditions  for  stability. 

We  conclude  this  section  v;ith  some  explicit  formulae 
for  the  various  terms  in  the  Taylor-expansion  of  the  energy 
expressions  for  our  examples  discussed  earlier.  In  the  case  of 
buckling  of  a  compressed  bar  we  have  from  (3.1.4) 
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P^[w]  =  0  , 


Pgtw]  =1  [J  Bw"2  -  I  Nw'2]dx  , 


p^Ew]  =  0  , 


Pj^Ew]  = 


Z 


E^  Bw'^w"^  -  ^  Nw'^jdx  , 


P3EW]  =  0  , 


PsEw]  = 


Z 


^  O 


Ei  Bw'\"^  -  ^  Nw'^Jdx  . 


(3.2.10) 

(3.2.11) 

(3.2.12) 

(3.2.13) 

(3.2.14) 

(3.2.15) 


For  the  stability  problem  of  a  three-dimensional 
elastic  body  under  dead  loading  we  have  from  (3.1.7)  and  (3.1.8) 

PlEu]  =  [  eJ  S^j(ui,/'ij,i)  -  x^u^]dv  -  I  p.u^dS  ,  (3.2.16) 

p 

PgCu]  =  1  (| 

+  '  (3.2.17) 

P3Eil]  =  I  ^ 

(3.2.18) 

P^Eu]  =  J  gE|  • 

(3.2.19) 


The  first  variation  of  the  energy  must  of  course 
vanish  In  both  examples^  because  the  fundamental  state  Is  a 
configuration  of  equilibrium.  In  our  first  example  this  fact 
is  immediately  confirmed  by  (3.2.10).  In  our  second  example 
we  have  to  keep  In  mind  that  the  distribution  of  initial 
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stresses  must  satisfy  the  equations  of  equilibrium 

Slj,J  Xi  =  0  (3.2.20) 

in  the  interior  of  the  body^  and  the  boundary  conditions 

SijHj  =  Pi  (3.2.21) 

on  the  surface  portion  S^.  The  vanishing  of  the  first  variation 
(3.2.16)  is  now  easily  confirmed  by  means  of  the  divergence 
theorem  (cf.  section  2.4  of  the  appendix)  and  the  geometric 
conditions  u^  =  0  on  S^.  Finally,  we  note  that  the  Taylor- 
expansion  of  the  energy  expression  (3.1.8)  for  the  second  example 
terminates  with  its  fourth  variation  (3. 2.19). 
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3 . 3  A  siifricient  condition  for  stability . 

It  will  first  be  shown,  by  means  of  a  simple  counter¬ 
example  of  a  v/ell-known  type  in  the  calculus  of  variations 
[e.g.  3,9],  that  a  positive  definite  second  variation  is  not 
always  a  sufficient  condition  for  a  proper  relative  minimijim. 

Our  counter-example  is  admittedly  somewhat  artificial  from  a 
physical  point  of  view.  Nevertheless  it  indicates  clearly  v;hy 
the  criterion  of  a  positive  definite  second  variation  may 
occasionally  fail  as  a  sufficient  condition.  It  also  points  the 
way  to  a  slightly  modified  and  uniformly  valid  sufficient  condi¬ 
tion  for  stability,  based  on  the  second  variation.  Our  modified 
criterion  reduces  to  the  condition  of  a  positive  definite  second 
variation  in  all  conceivable  practical  applications  (cf.  section 


In  our  counter-example  we  consider  an  expression  for 
the  energy  Increase  (3.1.1)  In  the  form  of  a  single  Integral 


P[u] 


n 


(x^u'^-u 


(3.3.1) 


o 

where  u(x)  is  the  single  displacement  component,  and  primes 
denote  differentiations  with  respect  to  x.  The  geometric  condi¬ 
tions  for  a  kinematically  admissible  displacement  field  are 
assumed  to  be  u(0)=u(l)=0. 

The  first  variation  of  (3.3.1)  is  identically  zero, 
and  the  second  variation 


rl 


x2u>2 


dx 


Pgtu]  = 


o 


(3.3.2) 
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is  evidently  positive  definite.  We  now  consider  a  one -parametric 
family  of  displacement  fields  u(x;e) 


0<x<e  <1  :u  =  e[x-2x^/e+x^/e'^]  , 


e  <  X  <  1 


:  u  =  0  . 


(3.3.3) 


Both  the  displacement  and  its  first  derivative  are  continuous 
and  tend  to  zero  for  e  — >  0.  The  energy  increase  (3.3.1)  for 
the  displacement  field  u(x;e)  is  easily  evaluated.  The  result 
is 

P[u(x;e)]  =  ^  e5  -  2_  ,  (3.3.4) 


where  the  first  term  is  due  to  the  second  variation.  This  result 
is  obviously  negative  for  all  e  in  the  range  0  <  e  <  1,  and  it 
follows  that  the  energy  increase ^  and  therefore  the  energy  itself ^ 
has  no  proper  minimum  for  u  =  0. 

The  failure  of  the  positive  definite  second  variation 
as  a  criterion  of  stability  in  the  present  example  is  of  course 
due  to  the  factor  x^  in  the  integrand  of  (3.3.2).  No  matter  how 
narrow  we  choose  the  neighbourhood  of  u  =  0^  that  is  no  matter 
how  small  we  choose  the  positive  bound  g*  for  the  modulus  of  u*, 
there  always  exist  displacement  fields  for  which  the  cubic  term 
in  (3.3.1)  is  the  dominant  term. 

An  obvious  way  out  of  our  difficulty  would  be  to 
require  that  not  only  the  second  variation  itself  but  also  its 
integrand  is  positive  definite.  Thus  it  might  be  ensured  that 
the  third-order  remainder  is  always  less  in  magnitude  than  the 
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second  variation  for  a  sufficiently  narrow  neighbourhood  of  the 
fundamental  state.  However^  our  mark  would  be  overshot  by  a 
wide  margin  by  this  more  sweeping  criterion.  Pew  problems  of 
elastic  stability  where  equilibrium  is  still  stable,  would  meet 
the  too  exacting  sufficient  condition  of  a  positive  definite 
integrand  of  the  second  variation  of  the  energy.  This  condition 
is  certainly  sufficient,  but  it  is  by  no  means  necessary. 

Our  purpose  is  actually  achieved  not  by  requiring  that 
the  second  variation  should  have  a  positive  definite  integrand, 
but  by  comparing  P2[hJ  with  an  auxiliary  homogeneous  quadratic 
fimctlonal  T^Lu]  whose  typical  Integrand  is  a  positive  definite 
quadratic  function  of  all  arguments  v/hich  appear  explicitly  In 
the  energy  Increase  (3.I.I). 

It  is  convenient  to  rewrite  the  typical  integrand  of 
the  second  variation  in  the  form  ' 

i  =  [|  ’  (3.3.5) 

where  the  constants  are  chosen  in  such  a  way  that  both 
quadratic  forms  in  (3.3.5)  are  positive  definite  in  all  arguments 
which  appear  explicitly  in  P(y;)  .  Without  loss  in  generality  we 
may  take  0^,^^  =  0  for  X  /  jx  .  The  non- vanishing  constants 
(X  not  summed)  are  then  all  positive.  The  sum  of  Integrals 
whose  Integrands  are  given  by  Is  denoted  as  the  auxiliary 

quadratic  functional  T2[u].  Since  both  terms  in  (3.3.5)  are 

in 

^The  basic  idea  of  splitting  up  the  typical  Integrand  in  the 
form  (3.3.5)  seems  to  be  due  to  Trefftz  [18],  although  he  does 
not  pursue  the  argument  to  its  end. 


/ 
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(3.3.6) 


positive,  we  have  the  inequality 

This  inequality  ensures  the  existence  of  a  greatest  lower  bound 
say  d,  for  the  left-hand  member  of  (3.3.6).  This 
greatest  lower  bound  is  characterized  by  two  properties  .  It 
ensures  the  inequality 


(3.3.7) 


for  every  non-vanishing  kinematically  admissible  displacement 
field  u.  On  the  other  hand,  it  implies  for  every  positive 
number  e,  the  existence  of  a  kinematically  admissible  displace¬ 
ment  field  v  satisfying  the  inequality 


Po[v] 


(3.3.8) 


We  shall  now  prove  that  a  positive  value  of  the 
greatest  lower  bound,  d  >  0,  is  a  sufficient  condition  for 
stability.  Prom  our  necessary  condition  (3«2.7),  and  from 

(3.3.8)  we  have  already  a  necessary  condition  for  stability  in 
terms  of  this  lower  bound,  viz.  d  ^  0.  In  the  case  of  a  semi- 
definite  second  variation  we  must  have  d=0.  We  know  from  our 
previous  discussion  that  the  additional  necessary  conditions 

(3.2.8)  and  (3.2.9)  have  to  be  satisfied  for  stability  in  that 
case.  It  follows  that  our  criterion  d  >  0  Is  the  sharpest 
possible  sufficient  condition  for  stability  in  terms  of  the 


second  variation  alone . 
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The  actual  proof  of  our  sufficient  condition  for 
stability  is  comparatively  simple.  From  (3.2.3)  and  (3.2.6)  we 
have 

P[u]  =  PgEu]  +  P^[u]  .  (3.3.9) 

Consulting  (3.2 .2) j  we  write  the  typical  Integrand  of  the  third- 
order  remainder  in  the  form 

4  =  ^4  ^  (3.3.10) 

where  is  an  abbreviation  of  the  expression  between  brackets . 
Remembering  our  particular  choice  for  the  integrands  of  T^Lu], 
we  may  write 

i  summed),  (3.3.11) 

Introducing  the  new  abbreviation 

— ——  -  c^  (X  not  summed),  (3.3.12) 

we  obtain  the  inequality 

l^yVl  i  •  (3.3.13) 

We  nov/  recall  that  the  third-order  derivatives  of  the  typical 
integrand  are  bounded.  Hence  we  may  make  the  quantities  in 
(3.3.10)  and  (3.3.13)  as  small  in  modulus  as  we  please  by  restrict¬ 
ing  our  attention  to  a  sufficiently  close  neighbourhood  of  the 
fundamental  state,  i.e.  by  choosing  the  positive  constants 
gjgSg'':»  ...  in  (3.1.2)  sufficiently  small.  It  follows  that  the 
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typical  Integrand  of  the  third-order  remainder  may  be  reduced  in 
absolute  value  to  an  arbitrarily  small  positive  fraction  of  the 
typical  integrand  of  T^Lu].  The  same  conclusion  then  applies 
a  fortiori  to  the  corresponding  integrals  or  sums  of  integrals. 

We  have  therefore  by  a  suitable  choice  of  the  neighbourhood 

1  ^  '  (3.3.1^) 

where  e  is  a  positive  nvmiber  as  small  as  we  please.  Prom  (3.3 • 7) ^ 
(3.3.9)  and  (3.3.1^)  we  have  finally 

P[u]  2  (d-e)T2[H]  .  (3.3.15) 

and  the  sufficiency  for  stability  of  the  criterion  d  >  0  has 
been  proven. 

It  will  be  observed  that  the  particular  choice  for  the 
typical  integrand  of  T^Tu]  is  immaterial,  provided  that  it  is 
positive  definite  in  all  arguments  appearing  in  the  energy 
increase  functional  (3.I.I).  We  might  even  relax  this  require¬ 
ment  slightly.  It  is  Indeed  sufficient  that  the  typical  inte¬ 
grand  of  1*2  [u]  is  positive  definite  in  all  argiaments  which  appear 
explicitly  in  the  third-order  remainder  (cf.  (3-3.10)).  The 
independence  of  our  criterion  of  the  particular  choice  of  typical 
integrand  for  T^Lu]  results  from  the  inequalities,  holding  for 
any  pair  of  positive  definite  quadratic  forms  [5,  ch.l] 


%yV  i  , 

(3.3.16) 

(3.3.17) 
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where  C  and  c  are  suitable  positive  numbers.  The  actual  value 
of  the  greatest  lower  bound  d,  defined  by  (3.3.7)  and  (3.3.8) ^ 
does  of  course  depend  on  the  choice  of  T^EuJ^  but  a  positive 
value  of  d  corresponding  to  T^Cu]  implies  a  positive  value  d* 

■K-  ‘K’ 

corresponding  to  T^Lu],  and  d=0  implies  d  =0,  and  vice  versa. 

We  shall  see  in  the  next  section  that  a  much  larger  freedom  in 
our  choice  of  T2[u]  even  exists  for  most  practical  applications. 

Finally,  it  may  be  worthwhile  to  return  briefly  to  our 
counter-example  (3.3.1)  concerning  the  criterion  of  a  positive 
definite  second  variation  of  the  energy  as  a  sufficient  condition 
for  stability.  It  is  easily  verified  that  our  modified  criterion 
based  on  the  second  variation  is  not  in  contradiction  to  the  non¬ 
existence  of  a  minimum  in  this  example .  If  we  take 

T  [u]  =  ui'^dx  ,  (3.3.18) 

o 

we  obtain  by  means  of  the  one-parametric  family  of  displacement 
fields  (3.3.3)  pi 

P„[u]  J  x^u'^dx 

d  =  g.i.b.  =  g.i.b.  -  <  lim  4  6“=^  =  0  .  (3.3.19) 

Pu'"^dx  -  e->o  ‘ 

Since  we  know  from  the  positive  definite  character  of  the  second 
variation  (3.3.2)  that  the  greatest  lower  bound  d  is  non-negative, 
it  must  necessarily  be  zero.  Our  modified  sufficient  condition 
for  stability  thus  fails  to  yield  a  decision,  and  no  contradiction 
arises  between  our  criterion  and  the  previously  established  non¬ 
existence  of  a  mlnimimi  of  (3.3.1)  for  u  =  0. 
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3.4  Reduction  to  a  variational  problem. 

Before  we  continue  our  discussion  it  may  be  helpful  to 
summarize  the  results  of  the  analysis  in  the  preceding  section. 

Let  TgLu]  denote  a  homogeneous  quadratic  functional  with  positive 
definite  integrands  in  all  arguments  which  appear  explicitly  in 
the  energy  increase  (3.I.I).  Let  d  denote  the  greatest  lower 
bound  of  the  quotient  of  the  second  variation  of  the  energy 
P2[u]  and  T^Lu] 

Pp[u] 

^  •  (3.4.1) 

A  necessary  condition  for  stability  of  our  fundamental  state  is 
then  d  >  0,  a  sufficient  condition  for  stability  is  d  >  0. 

Remarkably  simple  as  this  criterion  may  appear,  it  is 
very  difficult  to  apply  lonless  we  can  assume  that  the  greatest 
lov;er  bound  d  is  actually  attained  for  some  non-vanishing  kine¬ 
matically  admissible  displacement  field.  This  assxmiption  implies 
the  existence  of  a  mlnlmura  of  the  quotient  of  the  two  functionals 
^2^—^  and  T2[u].  We  emphasize  that  the  existence  of  such  a 
minimum  cannot  be  proved  in  complete  generality.  On  the  contrary, 
counter-examples  are  easily  constructed  in  v/hich  no  minimimi 
exists,  in  spite  of  the  existence  of  a  greatest  lower  bound.  A 
simple  example  of  this  type  is  provided  by  the  quotient  of  func¬ 
tionals  in  (3.3.19).  The  greatest  lower  bound  zero  in  this  case 
is  not  attained  for  any  non- vanishing  kinematically  admissible 

function  u(x) .  The  non-existence  of  a  minimum  is  again  due  to 
2 

the  factor  x  in  the  integrand  of  the  second  variation.  No 
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similar  difficulty  would  occur^  if  this  factor  were  replaced  by 
a  continuously  differentiable  function  of  x  with  a  positive 
minimum  in  the  closed  interval  0  _<  x  1. 

In  order  to  come  to  grips  with  our  stability  problem, 
we  shall  henceforward  assume  that  the  minimum  problem 

actually  has  a  solution,  specified  by  some  non-vanishing  kine¬ 
matically  admissible  displacement  field  u^ .  We  are  guided  in 
this  assumption  by  our  physical  intuition  chat  the  actual 
structure  of  the  Integrands  of  the  second  variation  p2[u]  in  a 
physically  meaningful  stability  problem  will  ensure  the  existence 
of  a  solution.  We  derive  some  support  for  our  intuition  from  the 
theory  of  small  vibrations  mentioned  below.  It  must  be  admitted, 
however,  that  a  complete  theory  would  require  a  verification,  at 
least  a  posteriori,  that  a  solution  of  our  minimum  problem  (3.^.2) 
actually  exists .  Such  an  existence  proof  is  indeed  available  for 
a  wide  class  of  problems  of  type  (3.^.2),  but  we  shall  not  pursue 
this  aspect  of  the  theory.  The  interested  reader  is  referred  to 
the  treatise  by  Courant  and  Hilbert  [6]  for  a  detailed  discussion 
of  existence  theory  in  variational  problems  of  quadratic  func¬ 
tionals  . 

V/e  may  now  rephrase  our  previous  criterion  of  stability 
in  terms  of  the  solution  of  our  minimum  problem  (3.‘*^-.2)  :  a 
necessary  condition  for  stability  of  our  fundamental  state  I  is 
^  sufficient  condition  for  stability  is  >  0. 
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decision  is  obtained  as  yet  In  the  critical  case  =  0,  which 
case  will  be  discussed  In  the  next  sections . 


An  interesting  consequence  of  our  assumption  of  the 
existence  of  a  solution  to  the  minimum  problem  (3.4.2)  is  that 
a  positive  definite  second  variation  of  the  energy  is  novj 
actually  a  sufficient  condition  for  stability.  For  if  a  solution 
to  problem  (3.4.2)  exists,  it  is  necessarily  positive  if  the 
second  variation  is  positive  definite.  Hence  we  return  to  the 
criterion  of  stability  in  its  customary  form,  which  we  had 
discredited  in  section  3.2.  A  casual  reader  may  nov/  well 
question  the  purpose  of  our  entire  discussion  in  section  3.3. 

We  underline  therefore  that  our  preceding  discussion  has  put  the 
conventional  criterion  of  a  positive  definite  second  variation 
as  a  sufficient  condition  for  stability  in  its  proper  perspective. 
It  is  justified,  if  and  only  if  a  solution  of  the  minimum 
problem  (3.4.2)  exists. 

A  second  important  consequence  of  our  assumption  of  a 
solution  to  (3.4.2)  is  that  we  may  now  replace  the  functional 


1*2 [u],  whose  integrand  is  positive  definite  in  all  argtiments 
appearing  explicitly  in  (3.1.1)^  by  any  positive  definite  homo- 


generous  quadratic  functional  T^[u],  and  consider  the  minimum 


problem 


=  Min. 


TgEu] 


(3.4.3) 


If  a  solution  to  (3.4.3)  also  exists,  it  must  obviously  have  the 
same  sign  as  the  solution  of  (3.4.2),  and  if  one  of  the  minima 
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or  is  zero,  the  other  one  must  also  be  zero.  The  corre- 

spending  displacement  fields  and  u^,  for  which  the  minima 
* 

and  are  attained,  are  of  course  not  the  same  in  general,  but 

if  ~  0,  we  may  also  Identify  the  displacement  fields  u^^ 

* 

and  u^ .  The  Increased  freedom  in  the  choice  of  an  auxiliary 
quadratic  functional  T^Lu]  is  convenient  for  many  applications. 
Whenever  we  use  this  greater  freedom,  we  shall  mark  this  fimc- 
tional  by  an  asterisk  as  in  (3.^.3).  freedom  is  of  course 

essentially  limited  by  the  requirement  that  the  functional 
T^luj  must  be  positive  definite,  even  if  its  typical  Integrand 
need  not  be  definite. 

As  a  first  Illustration  of  the  wider  class  of  function- 

* .  , 

als  T^Luj^  we  consider  the  expression  for  the  kinetic  energy  of 

our  mechanical  system.  Replacing  the  velocities  by  the  dlsplace- 

* 

ment  components  themselves,  we  obtain  a  functional  of  type  T^LuJ^ 
whose  typical  integrand  is  a  positive  definite  quadratic  function 
of  the  displacement  components.  The  minimum  property  (3.^.3) 
expresses  in  this  case  Rayleigh’s  principle  for  the  square  of  the 
fundamental  frequency  of  free  (small)  vibrations  of  our  system. 
The  displacement  field  for  which  the  minimum  is  attained,  is  the 
corresponding  fundamental  mode.  The  existence  of  such  a  funda¬ 
mental  mode  of  vibration  is  a  well-established  experimental  fact. 
It  ensures  the  existence  of  a  solution  to  the  minimum  problem 
(3. “^.3)  in  this  case,  and  it  supports  the  assumed  existence  of  a 
solution  to  our  original  minimum  problem  (3.^.2). 
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V/e  no\j  return  to  our  basic  minimum  problem  (3.^.2). 

Let  u^  again  denote  the  displacement  field  for  v/hich  the  minimum 
is  attained.  Let  C  denote  an  arbitrary  kinematically  admissible 
displacement  field.  Since  our  geometric  conditions  are  linear 
and  homogeneous,  the  displacement  field  also  kinematical¬ 

ly  admissible  for  any  value  of  the  constant  e.  The  minimum 
property  of  the  solution  of  (3.^.2)  is  then  expressed  by 


Tglu^+s^J  ^  ' 


(3.4.4) 


v/hence 

PgCui+eiJ  -  “iT2[u^+eC]  >  0  .  (3.4.5) 


V/e  introduce  a  convenient  notation  for  the  binomial 
expansion  of  any  functional  S^[u],  whose  typical  integrand  is  a 
homogeneous  polynomial  of  degree  m  in  the  components  of  the  dis¬ 
placement  u,  and  their  derivatives  u.  u.  etc.  We  write, 

if  u  =  ^w  , 


[  v+w  J 


m 

Z 

n=o 


(3.4.6) 


where  the  typical  integrand  of  S^^_^^^[v,w]  contains  all  terms 
in  the  binomial  expansion  of  degree  (m-n)  in  v  and  its  derivatives, 
and  of  degree  n  in  w  and  its  derivatives .  Accordingly  we  have 
the  identities 


=  S^(n,-n)  S^[v]  , 

(3.4.7) 


S„„[VjV/]  =  S  [w,v]  =  S  [v]  . 
mo  —  —  om  m  — 
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Applying  the  binomial  expansion  (3.^.6)  to  both  terms 
in  and  remembering  the  identity 

^  (3.^.8) 

we  obtain  from  (3. ‘^.5) 

e{P3_^[ui,^]  -  +  e^CPgE^]  -  t>>iT2[i]}  2  0  . 

(3.^.9) 

This  inequality  must  hold  for  arbitrary  values  of  the  constant  s . 
A  necessary  condition  for  (3.^.9)  to  hold  is  therefore 

P;^l[ui,C]  -  =  0  (3.4.10) 

for  every  kinematically  admissible  displacement  field  C.  In 
view  of  the  minimum  property  of  u)^  (3.4.2),  this  condition  is 
also  sufficient  to  satisfy  (3.4.9)j>  because  the  second  term  is 
obviously  always  non-negative. 

It  is  convenient  to  formulate  our  resulting  equation 
(3.4,10)  in  the  slightly  different  fom  that  u=co^  and  u=u^  are 
a  solution  of  the  equation 

=  0  ,  (3.4.11) 

holding  for  every  kinematically  admissible  displacement  field  Z. 
V/e  shall  call  this  more  general  equation  the  variational  equation 
associated  with  the  quotient  of  our  two  quadratic  functionals 
P^Lu]  and  T^Lu].  If  u  is  a  non-vanishing  solution  of  (3.4.11) 
for  some  value  of  the  parameter  u)  ,  we  obtain  from  the  equation 
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itself^ 


evaluated  for  ^  =  u, 

w  = 


Pgk] 

f^TuT  ’ 


O.'J.is) 


where  appropriate  use  has  been  made  of  (3.^.7).  An  even  more 
significant  property  of  a  solution  u  of  (3. ^.11)  is  that  it 
yields  a  stationary  value  of  the  quotient  in  (3-^. 12).  This 
statement  is  based  on  the  easily  verified  fact  that  the  quotient 
in  (3. ^.12)^  evaluated  for  any  kinematically  admissible  displace¬ 
ment  field  u+eC  ^  differs  from  (3. ^.12)  by  an  amount  which  is 
quadratic  in  e,  if  u  is  a  kinematically  admissible  non-vanishing 
solution  of  (3.‘^.ll)  .  The  absence  of  a  linear  term  in  0  is 
characteristic  of  the  stationary  property  of  cj  in  (3. ^.12). 

V/e  note  that  the  variational  equation  is  homogeneous 
and  linear  in  the  displacement  field  u.  Hence  it  will  always 
have  a  trivial  solution  u  =  0.  Non-vanishing  kinematically 
admissible  solutions  of  (3.'^.ll)  v/lll  be  called  (vector)  eigen- 
solutions  .  As  a  rule^  they  exist  only  for  special  values  of  the 
parameter  oj,  the  so-called  eigenvalues  which  coincide  with  the 
stationary  values  (3. ^.12).  In  view  of  the  hoinogeneity  of  eq. 

(3. ^.11)  any  eigensolution  must  contain  an  indeterminate  constant 
factor.  In  order  to  remove  this  ambiguity  we  shall  normalise 
our  eigensolutions  by  requiring  that  the  functional  T2[u]  takes 
a  fixed  positive  value.  It  is  usually  expedient  to  normalise  by 
talcing  T2[u]  equal  to  unity ,  and  we  shall  also  adopt  this  conven¬ 
tion.  The  normalised  eigensolutions  will  be  called  (vector) 
eigenfunctions .  We  shall  order  our  eigenvalues  in  an  increasing 
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sequence  ~  ^'^2  —  ^3  ~  *  ’  *  *  associated  eigenfunctions 

will  be  called  u^,  u^j  u^,  ...  This  nomenclature  ’  agrees  of 
course  with  our  previous  definition  of  and  u^,  although  v/e 
have  now  disposed  of  a  previously  indeterminate  constant  factor 
in  u ^ .  For  each  combination  of  eigenvalue  and  associated  eigen¬ 
function  we  have  the  identity 


An  important  property  of  the  eigenfunctions  u^  and  Uj^^ 
associated  with  different  eigenvalues  and  is  that  they  are 
orthogonal  with  reference  to  the  quadratic  functional  T2[u]. 
Orthogonality  of  two  kinematically  admissible  displacement  fields 
V  and  w  is  here  defined  by  the  relation 


Tiitv.w]  =  0  .  (3.4.14) 

We  note  that  the  positive  definite  character  of  T^Lu]  implies 
linear  Independence  of  the  displacement  fields  v  and  w  in  (3.4.14). 
The  proof  of  our  statement  with  respect  to  the  orthogonality  of 
u^  and  Uj^  follows  immediately  from  our  variational  equations 
(3.4.11)  j  satisfied  by  u^^  for  and  by  u^^  for  Taking 

^=Uj^  in  the  first  equation^  ^  =u^  in  the  second  equation^  and 
subtracting^  we  obtain  by  means  of  (3.4.7) 


]  =  0 


(3. ^.15) 


*1 

'The  reader  should  distinguish  carefully  between  our  present 
notation  for  the  eigenfunctions  u^^  (h=l,2j...),  each  of  ivhich 

represents  a  vector  field  with  local  displacement  components 
Uj^l^Uh2^ ^h3^  and  our  previous  notation  Uj(j=l,2j3)  for  the 

displacement  components  themselves. 
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Since  ^  Wj^,  v:e  have 


=  0  . 


(3.'^.l6) 


By  means  of  the  variational  equation  for  or  we  obtain  also 

Pll[Uh.Uk]  =  0  .  (3.4.17) 


Our  proof  of  orthogonality  of  tvzo  eigenfunctions 
and  u^  breaks  down,  if  the  associated  eigenvalues  are  equal , 

general  case  of  a  multiple  eigenvalue  is  described 

by 

“h  =  ‘“'(h+l)  =  •••  =  «(h+n)  •  (3. '^.18) 


The  equality  of  these  (n+l)  eigenvalues  implies  that  the 
variational  equation  has  (n+l)  linearly  Independent  elgensolutlons 
for  this  particular  value  of  w.  Let  — (h+l)  ^  *  *  *  •"— (h+n) 
the  associated  normalised  elgensolutlons  or  eigenfunctions.  Prom 
these  we  may  nov;  construct  a  new  set  of  (n+l)  eigenfunctions 
which  satisfy  our  orthogonality  condition  (3.4.16).  The  relation 
(3.4.17)  follows  then  again  immediately  from  the  variational 
equation . 


By  linear  combination  of  our  available  eigenfunctions 
j V(h^n)  construct  elgensolutlons  defined  by  the 
recurrent  scheme 


k-1 

^(h+k)  "  ^(h+k)  “kj  -(h+j) 


k=0.1,2,...n  (3.4.19) 


where  a,  ,  (k=l,2, . . .n; j  <  k)  are  constant  coefficients,  or  in 
more  explicit  form 
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Jih  =  :!!:h 

Ji(h+i)  =  ^(h+i)  ‘'lo  i4i 

-(h+2)  "  -(h+2)  °20  -h  “2I  -(h+l)  ’ 


ii(h+n)  =  ^(h+n)  +  “no  +  %(n-l)5i(h+n-l)  • 


From  the  construction  (3. ^.19)  it  is  immediately  obvious  that 
our  new  eigensolutions  are  again  linearly  independent.  We  deter¬ 
mine  the  coefficients  (k=l,2 ^ . .n; j  <  k)  recurrently  from  the 

requirement  that  each  eigensolution  shall  be  orthogonal 

to  all  its  predecessors  ,  j  <  k.  We  obtain  from  this 

condition 


a 


=  _  -(h+.l)  ^ 


kj  -  -  2T2lw(h+j) J 


(3. ^.20) 


Normalising  the  (n+l)  orthogonal  eigensolutions 
(j=0jl:,2,  ...,n)^  we  have  indeed  constructed  a  set  of  (n+l) 
mutually  orthogonal  elgenfimctlons  associated  with  the  multiple 
eigenvalue  (3.4.18). 


Henceforward  we  shall  assume  that  we  have  already 
orthogonalised  the  eigenfunctions  associated  with  a  multiple 
eigenvalue .  Orthogonality  is  then  a  general  property  of  any 
pair  of  eigenfunctions  which  enter  our  discussion. 

An  Important  consequence  of  the  orthogonality  of 
eigenfunctions  associated  with  different  eigenvalues  is  that 
all  eigenvalues  are  real.  This  theorem  is  easily  proved  in¬ 
directly.  If  it  were  false,  and  a  complex  eigenvalue  X=  +i|jL 
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would  exist,  its  associated  eigenfunction  would  also  possibly 
be  complex  u  =  y+iw.  Its  complex  conjugate  y-lw  would  then  also 
satisfy  the  variational  equation  (3.^.1l)  for  the  complex 
conjugate  X-ijx  of  u).  The  orthogonality  relation 

'^llUy+lw),  (y-lw)]  =  0  (3.4,21) 

is  incompatible  with  the  expansion  of  its  left-hand  member 

Tii[(y+lw)(y-iw)  ]  =  2T^[v]  +  >  0  .  (3.4.22) 

Our  contradiction  between  (3.4.21)  and  (3.4.22)  Implies  that  the 
assumption  of  existence  of  a  complex  eigenvalue  is  untenable. 

We  note  that  eigensolutions  of  the  variational  equation 
are  not  necessarily  real,  but  the  real  and  imaginary  parts  of  a 
complex  eigensolution,  taken  separately,  must  each  satisfy  this 
equation.  Any  complex  eigensolution  may  therefore  be  represented 
as  a  linear  combination  of  real  eigenfunctions  with  complex 
coefficients.  Without  loss  in  generality  we  may  restrict  our 
attention  to  the  real  eigenfunctions . 

A  further  important  consequence  of  our  orthogonality 
relations  is  that  we  may  write  any  kinematically  admissible  dis¬ 
placement  field  u  as  the  sirni  of  a  linear  combination  of  a  number 
of  eigenfunctions  and  a  residual  displacement  field  y  orthogonal 
to  these  eigenfunctions,  for  Instance 
k 

u  =  Z  a^Uj  H-  y,  where  T^^ [u^,y]  =  0,  j=l,2,...k.  (3.4.23) 
j  I 

In  order  to  prove  this  thoerem  we  write  the  required  orthogonality 
relations  for  v  in  the  form 
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(u  -  a.Uj]]  =  -  ^2^  ajT^^[Uj^,Uj  ]  = 


=  -  2aj^  =  0  ,  h  =  l,2,...k.  (3.'^.24) 

Our  purpose  is  achieved^  and  our  theorem  has  therefore  been 
proved,  if  we  take 

aj  =  I  ,  j=l,2,  ...k  .  (3.4,25) 

Except  for  the  first  eigenvalue  and  its  associated 
eigenfunction  u^,  our  ordered  sequence  of  eigenvalues  and  associ¬ 
ated  eigenfunctions  has  been  introduced  through  the  variational 
equation  (3.^.1l).  We  have  already  seen  that  each  eigenfunction 

u.  defines  a  stationary  value  w  of  the  quotient  (3. ^.12).  V/e 
J  j 

shall  nov7  prove  that  we  may  also  define  the  higher  eigenvalues 
and  their  associated  eigenfunctions  as  the  solutions  of  a  sequence 
of  minimum  problems  similar  to  (3.^*2),  The  existence  of  solutions 
to  the  minlmiAm  problems  in  this  sequence  is  of  course  tacitly 
assumed  in  this  alternative  definition.  The  only  modification 
in  our  sequence  of  minimum  problems,  in  comparison  with  our  basic 
problem  (3. ‘^.2),  consists  of  a  restriction  in  the  class  of 
competing  admissible  displacement  fields.  In  each  minimum 
problem  of  our  sequence  we  restrict  the  class  of  competing  kine¬ 
matically  admissible  displacement  fields  by  the  additional 
requirement  that  they  must  be  orthogonal  to  the  solutions  of  all 
preceding  minimum  problems.  Our  modified  definition  of  as 
the  solution  of  the  h-th  minimum  problem  reads  therefore 


3^ 


(3.4.26) 

subject  to  the  subsidiary  conditions  or  side  conditions 

Tii[Uj,u]  =  0  ,  j=l,2, ...(h-l)  .  (3.4.27) 

We  may  normalise  the  solution  u^.^  of  the  problem  by  requiring 

T2[uj^]  =  1  .  (3.4.28) 

Before  we  prove  that  the  minima  and  the  solutions 
u^  coincide  with  the  eigenvalues  and  associated  eigenfunctions 
as  defined  previously,  v;e  observe  that  the  sequence  of  minima 
(3. '4-. 26)  is  necessarily  non-decreasing 

^1  =  ^’^2  =  ^3  =  ‘  (3. ^.29) 

due  to  the  increasing  number  of  restrictions  imposed  on  the 
competing  displacement  fields.  We  also  remember  that  the  first 
minimum  (3. ^.26)  is  defined  v/lthout  side  conditions  and  is  there¬ 
fore  necessarily  equal  to  the  first  eigenvalue,  originally  defined 
by  (3.^.2).  Our  proof  proceeds  now  by  complete  induction.  If 
the  coincidence  of  (jj  and  u^  according  to  both  definition  is 
assumed  up  to  j=(h-l),  vie  shall  prove  that  it  also  holds  for 
j=h.  Since  we  are  assured  of  the  equivalence  for  J=l,  we  are 
in  possession  of  the  required  starting  point  for  a  proof  by 
Induction . 

Let  Uj^  denote  the  normalised  solution  of  (3. -^.26)  and 
(3. ^.27).  and  let  denote  an  arbitrary  kinematically  admissible 
displacement  field  subject  to  the  same  conditions  of  orthogonality 
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(3. ^.27),  i.e. 


=  0  >  J=l,2, . .  .(h-l)  .  (3^.30) 


The  minimum  property  of  (3. ^.26)  is  then  expressed  by  an  in- 
equalityj  similar  to  (3.^.^)> 


^TviJI  2 


for  arbitrary  values  of  the  constant  e .  We  rewrite  this  inequality- 
in  a  iorm  similar  to  (3.^.9) 


(3.4.32) 


Since  the  second  term  is  non-negative  for  all  j}  which  satisfy 
(3. ^.30),  a  necessary  and  sufficient  condition  resulting  from 
this  inequality  is 

-  ViiK'ii^  =  0  (3.4.33) 


for  every  kinematically  admissible  displacement  field  _r]  which 
satisfies  the  orthogonality  relations  (3.^.30). 

Our  equation  (3. ^.33)  bears  a  close  resemblance  to  the 
variational  equation  (3.^.1l),  but  it  is  clearly  distinct  from 
the  latter  equation  in  view  of  the  restrictions  imposed  on 
(3. ^.33).  We  shall  prove,  however,  that  (3. ^.33)  remains  valid, 
if  ri  is  replaced  by  any  kinematically  admissible  displacement 
field  C  not  subject  to  restrictions  (3.^.30).  It  will  be 
remembered  from  (3 -^.23)  that  we  may  always  write 
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h-l 

C  =  2  t.u.  +  j]  ,  (3.^.3^+) 

j=l  ^  ^ 

where  j}  satisfies  (3.^.30).  To  this  end  we  have  to  choose 
the  coefficients  tj  corresponding  to  (3. ^.25) 

tj  =  J  .  j=i,2,...(h-i).  (s.'f.as) 

We  now  replace  in  the  left-hand  member  of  (3. ^.33)  by  C  , 
substitute  from  (3.^.3^)  and  evaluate  the  resulting  expression. 
We  find 


(h-l)  ^ 


(3. ^.36) 


The  last  pair  of  terms  cancel  each  other  on  account  of  {3-^«33)* 
The  remaining  terms  vanish  separately.  This  is  obvious  for  the 
terms  T^^[uj^,Uj]  from  the  orthogonality  conditions  (3. ^.27) 
imposed  on  u^^.  We  also  know  that  the  variational  equation 
(3.^.1l)  is  satisfied  by  •  ( j=1^2,  .  . .  (h-l)  .  It  follows 

from  this  equation^,  by  taking  C  =  that  vanishes 

simultaneously  with  ] •  Hence  we  have  from  (3. ^.36)  and 

(3.4.33) 


-  ‘^h^ll^^h'^^  =  °  ’  (3.4.37) 


and  we  have  proved  that  and  u^,  defined  by  the  minimum  problem 
(3. ^.26)  and  (3.4.27),  are  solutions  of  the  variational  equation 
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This  completes  our  proof  of  the  equivalence  of  our 
two  alternative  definitions  of  and  as  eigenvalues  and 
eigenfunctions  of  (3.^.1!)  on  the  one  hand,  and  as  minima 
(3. ^.26)  and  minimising  displacement  fields  on  the  other  hand. 

Earlier  in  this  section  it  was  remarked  that  we  may 
replace  T2[u]  in  the  denominator  of  the  minimum  problem  (3.^.2) 

■it  r  1 

by  any  positive  definite  quadratic  functional  TgiuJ  v/hose  typical 
integrand  is  not  necessarily  definite.  An  Investigation  of  the 
modified  minimum  problem  (3. ■^.3)  yields  the  same  information  on 
the  stability  of  the  fundamental  state  as  an  investigation  of  the 
original  problem  (3.^.2).  The  analysis  of  the  latter  problem., 
starting  with  (3.^.^)  and  ending  with  equation  (3.^.3?)  may  be 
repeated  word  for  word  in  the  case  of  the  minimum  problem  (3.^.3) ^ 
if  we  replace  the  quadratic  and  bilinear  functionals  T2[u]  and 
T]^^[v,w],  the  eigenfunctions  u^^  and  the  eigenvalues  by  their 
corresponding  starred  counterparts .  If  we  take  in  particular  the 

*  r  1 

case  where  T^LuJ  is  the  denominator  in  Rayleigh’s  principle,  we 
recover  the  pertinent  well-knov/n  theorems  for  small  free  vibra¬ 
tions  about  the  configuration  of  equilibrium  in  state  I[e.g.5]. 

The  reader  will  have  noticed  that  our  discussion  in  the 
present  section  has  been  based  entirely  on  functional  relations, 
and  in  particular  on  the  variational  equation  (3.^i-.ll).  For  the 
basic  notions  involved  we  have  drawn  heavily  on  the  fundamental 
treatise  by  Courant  and  Hilbert  [5^6],  We  have  two  reasons  for 
our  preference  of  the  functional  approach  in  the  present  general 
discussion  over  the  more  familiar  eigenvalue  theory  of 
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homogeneous  boundary  problems  for  linear  differential  equations. 

Our  most  fundamental  reason  is  that  the  energy  criterion  for  the 
stability  of  equilibrium  of  continuous  bodies  represents  essential¬ 
ly  a  concept  of  functional  analysis.  Even  though  the  discussion 
of  the  variational  equation  (S.^.ll)  might  admittedly  be  replaced 
by  a  similar  discussion  of  the  equivalent  differential  equations 
and  boundary  conditions ^  the  continuation  of  the  investigation 
in  the  critical  case  oj^=0  in  sections  3.6  to  3.8  would  still 
require  a  return  to  the  functional  formulation.  A  second  reason 
for  our  preference  is  the  far  greater  simplicity  of  the  discussion 
in  terms  of  functional  analysis,  in  comparison  with  a  discussion 
of  the  equivalent  differential  equations  and  boundary  conditions. 
All  relevant  conditions  are  embraced  by  a  single  variational 
equation.  Once  the  fev;  required  basic  concepts  of  functional 
analysis  have  been  grasped,  the  analysis  is  conspicuously  brief. 

Even  if  our  purpose  of  a  general  discussion  is  served 
best  by  a  formulation  in  terms  of  concepts  of  functional  analysis, 
a  better  understanding  of  its  implications  may  well  be  furthered 
by  a  transcription  into  terms  of  the  more  familiar  concepts  of 
differential  equations  and  boundary  conditions .  V/e  prefer  to 
postpone  such  a  transcription  until  the  end  of  the  next  section, 
where  we  shall  derive  from  our  variational  equation  (3.^.1l)  the 
differential  equations  and  boundary  conditions  for  neutral 
equilibrium  for  the  tv/o  examples  discussed  previously. 
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3.5  The  stability  limits  neutral  equilibrium  and  buckling  modes. 

In  the  last  section  we  have  seen  that  the  conditions 
for  the  stability  of  equilibrium  in  the  fundamental  state  I  may 
be  formulated  in  terms  of  the  solution  of  the  minimum  problem 
(3.^.2).  A  necessary  condition  for  stability  is  >  0^,  a 
sufficient  condition  is  >  0. 

No  decision  has  as  yet  been  obtained  in  the  limiting 
case  =  0.  This  case  will  be  called  a  critical  case  of  neutral 
equilibrium,  and  the  equilibrium  is  said  to  be  at  the  stability 
limit  in  this  case.  The  stability  limit  is  characterized  by 


CJ 


1 


^an. 


0  . 


(3.5.1) 


The  normalised  minimising  displacement  field  or  associated 
eigenfunction  u^  will  be  called  the  critical  buckling  mode.  We 
may  also  describe  the  case  of  equilibrium  at  the  stability  limit 
in  terms  of  the  semi -definite  second  variation  alone 


PgCu]  >  0 ,  °  ’  (3.5.2) 

where  u^  is  the  critical  buckling  mode.  Finally,  our  general 
variational  equation  (3.4.11)  takes  the  simplified  form  that  the 
buckling  mode  u^  is  a  solution  of  the  variational  equation 

=  0  ,  (3.5.3) 

holding  for  every  kinematically  admissible  displacement  field  C. 

In  the  continued  investigation  of  stability  in  a 
critical  case  of  neutral  equilibrium,  to  be  discussed  in  sections 
3.6  to  3.8,  it  is  essential  to  know  first  all  possible  critical 
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buckling  modes ,  l.e.  all  displacement  fields  u  for  which  the 
second  variation  ^  vanishes.  In  order  to  answer  the  question 
whether  the  second  variation  may  also  be  zero  for  a  displacement 
field  u  which  is  not  a  multiple  of  u^^  v/e  write  (cf .  (3.4.23)) 

u  =  au^+v  ^  where  0  .  (3*5.^) 

We  have  now  from  (3.5.2)  and  (3.5.3) 

^2^-^  =  ""  a^P^Eu^]  +  ap^^[u^,v]  +  P^tv]  =  P2[v]  . 

(3.5.5) 


Hence,  if  P2[liJ  vanishes  for  a  displacement  field  u  which  is  not 
a  multiple  of  the  critical  buckling  mode  u^,  it  must  also  vanish 
for  the  displacement  field  v  orthogonal  to  this  buckling  mode, 
defined  by  (3.5.4). 

We  now  recall  our  sequence  of  minimum  problems  (3.4.26), 
(3.4.27) ^  in  particular 


(*) 


Min. 


(3.5.6) 


under  the  side  condition 


•  (3.5.7) 

If  is  positive,  no  displacement  field  v  orthogonal  to  u^  exists 
for  which  the  second  variation  is  zero .  In  this  case  the  critical 
buckling  mode  u^  is  called  a  simple  critical  buckling  mode.  On 
the  other  hand,  tx)2=W;j^=0^  a  second  critical  buckling  mode  exists 
in  the  form  of  the  eigenfunction  associated  with  the  second 
eigenvalue  ^^=0. 
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In  the  latter  case  of  a  so-called  multiple  critical 
buckling  mode,  we  must  of  course  continue  our  Investigation.  We 
now  have  to  answer  the  question  whether  displacement  fields  u 
exist  which  are  linearly  independent  of  both  critical  buckling 
modes  and  for  v:hlch  the  second  variation  again  vanishes .  By  a 
similar  argument  as  before  we  may  prove  that  the  answer  is 
provided  by  the  solution  of  the  third  minimimi  problem  in  the 
sequence  (3.4.26)^  (3.4.27).  If  is  positive,  no  displacement 
field  u  exists  which  is  linearly  independent  of  u-,  and  u^,  and 
for  v/hlch  P2[h^  zero.  Evidently  v;e  have  to  continue  our 
investigation  again,  if  ^^=0. 

Our  results  may  be  summarized  by  a  general  prescription 

of  the  procedure  to  be  followed  in  the  analysis .  The  sequence  of 

minimum  problems  (3.4,26),  (3.4.27)  has  to  be  investigated  until 

a  positive  minimum,  say  Is  obtained,  or  at  least  has  been 

shown  to  exist.  The  most  general  displacement  field  u  for  which 

the  second  variation  12^—^  vanishes  is  then  given  by  a  linear 

combination  of  h  critical  buckling  modes 

h 

u  =  S  a.u.  .  (3.5.8) 

j=l  ^  ^ 

It  will  be  convenient  to  refer  to  a  coefficient  a.  in  (3.5.8)  as 

3 

the  amplitude  of  the  buckling  mode  u^ . 

It  v/lll  be  noticed  that  the  vanishing  of  the  second 
variation  P2^— ^  ^  displacement  field  (3.5.8)  in  the  case  of 

a  multiple  eigenvalue  ...  =o^=0,  is  an  immediate  consequence 

of  the  discussion  of  the  variational  equation  (3.4.11)  in  the 
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previous  section.  It  had  not  been  proved  previously,  however, 
that  the  second  variation  cannot  vanish  for  other  displacement 
fields .  The  additional  analysis  in  the  present  section  was  Indeed 
required  to  establish  this  important  feature. 

The  conventional  approach  to  the  theory  of  elastic 
stability  starts  from  the  concept  of  neutral  equilibrium 
[e.g.  1,2,15].  In  this  approach,  a  fundamental  state  of  equi¬ 
librium  is  called  neutral,  if  additional  configurations  of  equi¬ 
librium  exist  which  are  obtained  from  the  fundamental  state  by 
means  of  infinitesimal  additional  displacements.  The  differential 
equations  and  boundary  conditions  which  express  the  existence  of 
such  adjacent  configurations  of  equilibrium,  are  linear  and 
homogeneous  in  the  additional  displacements  from  the  fundamental 
state .  These  equations  are  the  complete  equivalent  of  our 
variational  equation  (3.5»3)  .  This  statement  is  easily 

verified  by  means  of  (3.2.3)  and  (3.2.6).  Neglecting  higher 
order  terms  in  the  potential  energy  in  an  adjacent  state  II, 
obtained  by  the  infinitesimal  displacement  field  u  from  the 
fundamental  state  I,  we  have  for  the  energy  in  state  II 

Pii  =  Pj  +  Patu]  (3.5.9) 

The  condition  of  equilibrium  in  this  state  is  given  by  the  condi¬ 
tion  for  a  stationary  value  of  the  energy,  i.e.  the  variational 
equation  (3.5.3)  holding  for  any  kinematically  admissible  dis¬ 
placement  field  C . 

in - -  . 

‘  In  spite  of  occasional  claims  to  the  contrary  [e.g.  13 ]>  the 
approach  through  the  concept  of  neutral  equilibrium  is  subject 
to  the  same  restrictions  as  the  energy  criterion :  a  Justifica¬ 
tion  of  the  neutral  equilibrium  approach  to  stability  theory 
is  possible  only  for  conservative  systems  (cf.  sections  1.4 
and  2  .  .  . . )  . 
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The  equation  of  neutral  equilihrlum  (3.5.3)  is 
satisfied  if  one  of  the  eigenvalues  of  our  variational  equation 
(3.4.11)  is  zero.  It  is  immaterial  from  the  viewpoint  of  neutral 
equilibrium  according  to  its  classical  definition^  whether  or  not 
the  eigenvalue  zero  is  the  smallest  eigenvalue  of  (3.4.11).  In 
the  first-mentioned  case  we  have  and  vie  are  dealing  with 

a  critical  case  of  neutral  equilibrium  as  defined  earlier  in  this 
section.  In  all  other  cases  of  neutra.!  equilibrium  we  have 
<  0^  and  we  know  already  from  our  energy  considerations  on 
the  basis  of  the  second  variation  alone  that  equilibrium  is 
actually  unstable.  It  is  somewhat  unfortunate  that  the  conven¬ 
tional  notion  of  neutral  equilibrium  does  not  distinguish  clearly 
between  the  critical  case^  where  the  second  variation  of  the 
energy  alone  is  Incapable  of  supplying  a  decision  on  stability, 
and  the  other  cases,  where  the  second  variation  has  already 
established  instability.  The  critical  case  of  neutral  equilib¬ 
rium  is  by  far  the  most  important  one,  and  whenever  there  is  no 
danger  of  confusion  with  other  cases  of  neutral  equilibrium,  we 
shall  emit  the  words  "critical  case"  in  referring  to  the  neutral 
equilibrium  under  consideration. 

V/e  now  return  one  again  to  the  two  examples  which  have 
been  discussed  before  as  Illustrations  of  the  general  theory. 

V/e  shall  derive  for  these  examples  the  differential  equations  and 
boundary  conditions  which  ai’e  the  equivalent  of  the  variational 
equation  of  neutral  equilibrium  (3.5.3). 
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In  the  case  of  the  bar  under  compressive  end  loads, 
the  second  variation  of  the  energy  is  given  by  (3.2.11).  The 
geometric  boundary  conditions  are  in  this  case  w(0)=w(>^)=0 .  The 
explicit  fom  of  variational  equation  (3.5.3)  is 

'i 

Pll[w,C]  =  [Bw"C”-Nw'C']cax  =  0  ,  (3.5.10) 

o 

holding  for  every  tv;ice  continuously  differentiable  function 
C(x)  which  vanishes  for  x=0  and  x=^ .  We  proceed  to  remove  the 
derivatives  of  the  arbitrary  function  C(x)  under  the  sign  of 
integration.  This  purpose  may  be  accomplished  through  integration 
by  parts,  if  it  is  assumed  that  the  solution  of  (3.5*10)  has 
continuous  derivatives  up  to  and  including  the  fourth  order. 

The  result  of  a  repeated  integration  by  parts  is  then 

Bw"C'  -  [(Bw")'  +  Nw']C  +  [(Bw")"  +  (N\ir')']C  dx  =  0  , 

o  0  o 

(3.5.11) 

The  first  txvo  terms  in  (3.5.11) j  which  depend  only  on  the 
boundary  values  of  the  functions  concerned  and  their  derivatives 
will  be  called  boundary  terms.  In  view  of  the  kinematic  condi¬ 
tions  imposed  on  C(x),  that  is  C(0)=C(>^)=0,  the  second  boundary 
tem  vanishes  identically.  Since  both  N  and  B  are  constants, 
equation  (3.5.11)  reduces  to 

Bw”C'  +  J  [Bw""  +  Nw"]C  dx  =  0  .  (3.5.12) 

o  o 

Tlie  remaining  argument  is  based  on  a  fundamental  lemma 
of  the  calculus  of  variations  [5^ch.4].  Ilie  importance  of  this 
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lemma  for  our  subsequent  analysis  v/arrants  a  comprehensive  state¬ 
ment  and  proof  here.  Let  x  denote  a  set  of  n  independent 
variables  x^(i=lj2, . . .n) ^  and  let  f(x)  denote  a  function  of  these 
variables,  continuous  in  a  closed  domain  D  of  n-dlmensional 
x-space.  Let  (k=0,l,2, ..)  denote  the  class  of  functions  C(x) 
in  the  closed  domain  D  which  have  continuous  partial  derivatives 
up  to  and  including  the  order  k.  This  class  may  possibly  also 
be  restricted  by  certain  homogeneous  boundary  conditions  ^  on  the 
boundary  of  D.  If  the  relation 

.  f(x)i:(x)dx  =  0  (3.5.13) 

D 

holds  for  every  function  C(x)  of  class  C^,  it  follows  that 

f(x)  =  0  (3.5.14) 


everywhere  in  D.  An  indirect  proof  of  this  basic  lemma  is 
easily  given.  Suppose  it  were  false,  f(x)  would  be  non-zero, 
say  positive,  at  some  point  y  in  D.  The  continuity  of  f(x) 
implies  that  this  function  is  then  also  positive  in  some  n- 


dlmensional  sphere  of  radius  r  about  the  center  V/e  may  choose 
a  function  class  defined  by 


i:(x) 

ax) 


[r^  -  Z  (x  -y for  Z  (x  -y  )2  < 

i=l  1  i  i=l  3.  1  _  I 

r  (3.5.1^ 

0  everyiA.^here  else  in  D. 


The  theorem  remains  valid,  if  non-homogeneous  boundary  condi¬ 
tions  are  imposed  on  ^  complete  proof  is  more 

complicated  in  this  case. 
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This  function  C(2i)  positive  inside  the  sphere  of  radius  r 
about  zero  everywhere  else.  The  integrand  of  (3.5.13) 

is  therefore  positive  inside  this  sphere,  and  zero  everywhere 
else  in  D.  The  Integral  in  the  left-hand  member  of  (3.5.13)  is 
now  also  positive,  and  we  have  arrived  at  a  contradiction.  It 
follows  that  our  assiomption  of  a  non-zero  value  of  f(x)  SLt  the 
point  y;  is  untenable,  and  the  lemma  has  thus  been  proven. 

In  order  to  apply  the  fundamental  lemma  to  (3.5.12), 
we  consider  the  class  of  twice  continuously  differentiable  func¬ 
tions  C(x)  for  which  both  the  function  and  its  first  derivative 
vanish  at  x=0  and  x=^.  The  boundary  term  in  (3.5.12)  then  dis¬ 
appears,  and  we  obtain  the  well-known  homogeneous  and  linear 
differential  equation  of  neutral  equilibrium  for  a  bar  under 
compressive  end  loads 

Bw""  +  Nw”  =  0  .  (3.5.16) 

The  integral  in  (3.5.12)  now  vanishes  identically.  Removing  our 
artificial  restriction  C*(C)  =  C'(^)  =  0,  which  is  not  imposed 
by  any  geometric  condition  of  the  problem,  we  infer  that  the 
cofactor  of  C’  1^  the  boundary  term  must  also  vanish.  Hence  we 
obtain  the  homogeneous  and  linear  dynamic  or  natural  boundary 
conditions 

x=0  and  x=^ :  Bw"  =  0  ,  (3.5.17) 

in  addition  to  the  homogeneous  and  linear  geometric  boundary 
conditions 


x=0  and  x=i :  w  =  0 


3 


(3.5.18) 
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imposed  on  all  competing  functions  w(x) . 

The  solution  of  the  equation  and  boundary  conditions 
for  neutral  equllibriimi  is  extremely  simple  in  the  present 
example  because  (3.5.16)  is  a  homogeneous  linear  differential 
equation  with  constant  coefficients .  Non-trivlal  solutions 
which  satisfy  the  boundary  conditions  exist  only  for  special 
values  of  the  compressive  load,  the  so-called  buckling  loads 
[17,  i2.2] 

N  =  N  ,  k=l,2,...  (3.5.19) 

K  jfd. 


The  associated  solutions  are  called  buckling  modes .  They  are 
given  by 

w  =  sin  kTLx/>^  ,  (3.5.20) 

where  is  an  arbitrary  constant.  In  order  to  decide  whether 
the  neutral  equilibrium  at  a  buckling  load  13.5.19)  is  a  critical 
case  of  neutral  equilibrium,  we  have  to  investigate  the  second 
variation  (3.2.11)  itself.  This  is  again  quite  simple  in  the 
present  example.  Since  the  competing  deflections  w(x)  in  (3.2.11) 
are  at  least  twice  continuously  differentiable,  the  Fourier  sine 
series  of  an  arbitrary  deflection  function  satisfying  (3.5.17) 

and  (3.5.18) 

00 

w(x)  =  S  V/,  sin  kTix/,^  (3.5.21) 

k=l  ^ 


may  be  differentiated  twice  term  by  term.  The  second  variation 
(3.2.11)  may  then  be  written  in  the  form 


2  00  o 

=  77^  E  k  [ - ^  -  N]  W! 

^  k=l  $r 


k 


(3.5.22) 
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It  is  positive  definite  for  N  <  serai -definite  for  N=N^  and 
indefinite  for  N  >  .  Hence  N=N^  defines  a  critical  case  of 

neutral  equilibrium^  and  this  load  is  the  critical  buckling  load 
or  Euler  buckling  load.  At  the  higher  buckling  loads  N^(k  >  2) 
we  have  unstable  cases  of  neutral  equilibrium. 

We  now  turn  to  our  second  example,  the  three-dimensional 
elastic  body  under  dead  loads .  The  second  variation  of  the 
potential  energy  in  this  case  is  defined  by  (3.2.17)^  and  the 
geometric  boundary  conditions  for  a  kinematically  admissible 
displacement  field  are  u=0  on  S^.  It  is  convenient  to  Introduce 
the  concept  of  a  (symmetric)  fictitious  linear  strain  tensor 


'’ij  •  • 


(3.5.23) 


This  fictitious  strain  tensor  should  be  clearly  distinguished 
from  the  actual  strain  tensor  defined  by  (3.1.?) >  which 

describes  the  actual  deformation  of  the  body  in  the  transition 
from  state  I  to  state  II.  The  fictitious  strain  tensor  coincides 
with  the  actual  strain  tensor,  if  and  only  if  the  displacement 
field  u  from  the  fundamental  state  I  to  state  II  is  Infinitesimal . 
Likewise,  it  is  convenient  to  denote  the  second  term  in  the 
Integrand  of  (3.2.17) 


A  =  G[0,  ,0,  ,  + 


r(eKJ‘'] 


(3.5.2Jf) 


Ij  Ij  l-2v^"hh' 
as  the  fictitious  elastic  strain  energy  density,  associated  with 
the  fictitious  linear  strain  tensor  through  Hookers  law.  It 
coincides  with  the  actual  elastic  strain  energy  density,  if  and 
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cnly  if  the  fundamental  state  I  is  the  undeformed  state  and  the 
displacement  field  u  from  this  fundamental  state  to  state  II  is 
infinitesimal.  Finally^  v/e  introduce  the  (symmetric)  fictitious 
stress  tensor,  associated  with  the  fictitious  linear  strain 
tensor  through  Hooke's  law. 


= 


(3.5.25) 


where  is  the  unit  tensor  (gj^j=0  for  ®ij”^  • 

Here  again  the  fictitious  stress  tensor  coincides  v;ith  the 
actual  stress  tensor,  if  and  only  if  the  fundamental  state  I  is 
the  undeformed  state  and  the  displacement  field  u  from  the 
fundamental  state  to  state  II  is  infinitesimal . 

After  these  preliminaries  we  may  write  down  the 
explicit  form  of  our  variational  equation  (3.5*3)  for  neutral 
equilibrium 

^  °  ’  (3.5.26) 


where  appropriate  use  has  been  made  of  the  symmetry  of  both  the 
initial  stress  tensor  3^^  and  the  fictitious  stress  tensor 
Equation  (3.5.26)  must  hold  for  every  once  continuously  differ¬ 
entiable  kinematically  admissible  displacement  field  C.  We  may 
again  remove  the  derivatives  of  the  components  from  our 
integral  by  applying  Green’s  theorem.  To  this  end  we  have  to 
assume  that  the  displacement  fields  u  which  satisfy  (3*5.26)  are 
twice  continuously  differentiable.  On  this  assumption  v/e  obtain 
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j 


=  0 


(3.5.27) 


where  n  is  the  unit  outward  normal  vector  to  S.  Restricting  our 
attention  first  to  displacement  fields  C  which  vanish  on  Sy  we 
obtain  from  the  fundamental  lemma  of  the  calculus  of  variations 
the  homogeneous  and  linear  differential  equations  of  neutral 
equilibrium  [19] 


+  "hjtj  =  ° 


(3.5.28) 


Returning  to  (3.5.27)  and  remembering  that  no  geometric  boundary 
conditions  are  imposed  on  ^  on  the  part  S  of  the  surface ^  we 
obtain  the  homogeneous  and  linear  dynamic  or  natural  boundaiT 
conditions 


^  °  •  (3.5.29) 


In  addition  to  these  we  have  the  geometric  boundary  conditions 

u^  =  0  on  .  (3.5.30) 


A  slightly  different  form  of  the  differential  equations 
(3.5.28)  may  be  obtained  by  appropriate  use  of  the  original  equa¬ 
tions  of  equilibrium  (3.2.20)  for  the  initial  stresses  in  state 
I.  The  modified  equations  read 


d,  .  .  H-  S,  .u,  .  .  -  X.u,  .  =0 

bj, j  ij  h,ij  i  h,i 


(3.5.31) 


The  boundary  conditions  (3.5.29)  may  be  reduced  in  a  similar 

hy  inenns  of  the  oidgiiial  boni)dai*y  conditions  (3.2.21) 

for  S.  .  on  S  ^  to 
1  j  p 

PA,i  = 


(3.5.32) 
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The  dorl vat-ion  of  the  equations  and  boundary  conditions 
for  neutral  equilibrium  suffices  for  our  present  purpose  to 
demonstrate  the  appropriate  manipulation  of  the  general  varia¬ 
tional  equation  (3 *5 *3)  in  the  example  of  an  elastic  body.  The 
reader  is  referred  to  chapter  6  for  a  more  detailed  discussion 
of  neutral  equilibrium  in  this  case. 
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3 . 6  Stability  In  a  critical  case  of  neutral  equlllbrlijm . 

For  the  further  Investigation  of  stability  in  a  critical 
case  of  neutral  equilibrium  it  will  be  convenient  to  consider 
first  the  case  of  a  simple  critical  buckling  mode  u^ .  The  solu¬ 
tion  minimum  problem  (3.5 *6),  (3.5.?)  is  then  positive . 

Writing  an  arbitrary  kinematically  admissible  displacement  field 
u  in  the  form  (3. ^.23) 

u  =  au^+v  j  where  ^  ,  (3.6.1) 

we  are  thus  ensured  of  the  inequality 

^2^^^  2  '  (3.6.2) 

Where  to^  is  a  positive  constant. 

Our  reason  for  v/ritlng  an  arbitrary  kinematically 
admissible  displacement  field  in  the  form  (3.6.1)  is  our  conjec¬ 
ture  that  only  ’’small"  deviations  v  from  a  multiple  au^  of  the 
buckling  mode  need  to  be  considered.  In  fact,  it  is  not  difficult 
to  prove  by  an  argument  similar  to  our  discussion  in  section  3.3 
that  (3.6.2)  ensures  a  positive  energy  increment  (3.I.I)  for 
sufficiently  small  values  of  the  amplitude  a  in  (3.6.I),  provided 
that  a”  T^[v]  has  a  non-vanishing  limit  for  a  — >  0.  Hence  a 
further  investigation  is  only  required  for  those  deviations  v 

_P 

from  a  multiple  au^  of  the  buckling  mode  for  which  oT 
approaches  zero  if  a  tends  to  zero.  A  detailed  proof  of  these 
statements  may  be  omitted  here  because  the  following  analysis  does 
not  presuppose  any  restriction  on  the  order  of  magnitude  of  v  in 


comparison  with  au^ . 
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In  section  3»1  the  existence  of  some  positive  constants 
with  the  property  that  the  inequalities  (3.1.2) 
ensure  a  non-negative  energy  increase  (3.I.I)  for  every  kine¬ 
matically  admissible  displacement  field  u,  was  established  as  a 
necessary  and  sufficient  condition  for  the  stability  of  equilib¬ 
rium  in  the  fundamental  state  I.  A  slightly  modified  version  of 
this  criterion  is  appropriate  for  displacement  fields  written  in 
the  form  (3.6.I).  A  necessary  and  sufficient  condition  for 
stability  is  the  existence  of  some  positive  constants  A,g,g ' ^g”, . . . 
with  the  property  that  the  inequalities 

I  a|  1  A  ,  I  v^l  1  g  ,  I  v^^  ^  g*  ,  I  1  g”  ,  ...  (3.6.3) 

ensure  a  non-negative  energy  increase  (3.I.I) 

P[au^+y]  2  0  .  (3.6.4) 

The  equivalence  of  this  alternative  form  with  the  original 
criterion  is  easily  proved.  The  inequalities  (3.6.3)  evidently 
imply  the  inequalities  (3.1.2),  and  it  only  remains  to  show  that 
the  converse  is  also  true.  From  (3.4.25)  we  have 

a  =  -I  ,  (3.6.5) 

and  it  follows  that  a  is  finite  since  u  and  its  derivatives  are 
bounded.  Returning  to  (3.6.I),  we  see  that  v  and  its  derivatives 
are  then  also  bounded  simultaneously  v;ith  u  and  its  derivatives. 
Hence  the  inequalities  (3.6.3)  are  Indeed  equivalent  to  the 
inequalities  (3.1.2). 


The  arguments  of  the  typical  integrand  P(y;)  of  P[u], 
defined  in  section  3.2^  may  be  written  in  a  form  similar  to 
(3.6.1) 

y^  =  ay^  +  ,  or  ^  =  a^^^  +  z  ,  (3.6.6) 

where  y^  represents  these  argments  corresponding  to  the  critical 
buckling  mode  u^ .  We  apply  Taylor*s  formula  in  order  to  obtain 
an  expainsion  of  the  typical  integrand  about  the  set  of  arguments 
•  Teminating  this  expansion  at  the  third  order  remainder^ 
we  have 

+  §  ,  (3.6.7) 

where  the  arguments  in  the  third  order  derivative  are  given  by 
ay^  +  =  ay^  +  (X  not  summed;  0  ^  0^  <  l)  .  (3.6.8) 

We  identify  the  first  term  in  (3.6.7)  as  the  typical  integrand  of 
P[au^].  Collecting  the  integrals  of  the  first  terms  in  the 
Integrands  (3.6.7),  and  remembering  (3.2.3) ^  (3.2.6)  and  (3.5.2), 
we  may  express  this  fact  by  the  symbolism 

=»  +  a^P2^[u^]  +  0(a5)  .  (3.6.9) 

Next  we  expand  the  first  and  second  derivatives,  which 
appear  in  (3.6.7)  at  the  set  of  arguments  ay^  ,  by  means  of 
Taylor *s  formula 
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Xjx 


+ 


)  X 


(3.6,10) 


F 


,x,<“ii>  -  ',x„<”>  *  ■  <0  i  =  i  I)-, 


.  X[JL 


Xjj.'' 


(3.6.11) 

V/e  collect  the  integrals  whose  Integrands  are  linear  in  z  and 
homogeneous  of  degrees  zero  to  four  in  the  amplitude  a  of  the 
buckling  mode.  Remembering  (3.2.6),  (3.'*^-.6)  and  (3.5.1)  ^  we 
note  the  symbolic  relations 


zV^x(O) 


P^[v]  = 


M  ^''"Myyiy>,x,xvpd' 


(3.6,12) 

:]  =  0  , 

(3.6,13) 

(3.6.14) 

>3i[Hi.v]  , 

(3.6.15) 

(3.6.16) 

Finally,  we  identify  the  term  which  is  homogeneous  and  quadratic 
in  z,  and  independent  of  a  as  the  typical  integrand  of  ^2^— 
i.e.  the  symbolic  relation 

i  ^^^(0)  V^[v]  .  (3.6.17) 

All  terms  in  (3.6.7)  have  now  been  identified,  except 
the  last,  cubic  term  in  z  and  that  part  of  the  second,  quadratic 
term  in  z  whose  coefficient  is  given  by  the  last  term  in  (3.6.11). 
We  write  these  terms  together  in  the  form  A^^z^z^^,  where  the 
coefficients  A^^^^  are  specified  by 
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\  =  i  ^  i  •  (3-6-18) 

Denoting  the  sum  of  integrals  of  these  remaining  terms  by  R2[v]> 
we  have  the  symbolic  relation 

Rgtv]  .  (3.6.19) 

By  means  of  (3.6.9),  (3.6.12)  -  (3.6.1?)  and  (3.6.19) 
we  may  now  write  for  the  energy  Increase  (3.1.1)  in  the  kine¬ 
matically  admissible  displacement  field  (3.6.1) 

P[aui+v]  =  a3p2[u^]  +  a^P|^[u^]  +  a^Pg^^Eu^.v]  +  a^P^^Eu^^v]  + 


==^  RoEv]  . 


+  a^Pi^^Eu^jv]  +  PgEv]  +  RgEv]  +  O(a^)  . 


(3.6.20) 


If  we  restrict  our  displacement  field  to  a  multiple  of 
the  buckling  mode^  i.e.  if  v;e  take  v  =  0^  we  obtain 


P[au^]  =  a^P^[u^]  -i*  a'^PiiLu^]  +  O(a^)  .  (3.6.21) 


5' 


In  view  of  the  condition  of  a  non-negative  energy  increase  in  a 
stable  fundamental  state,  we  recover  the  necessary  conditions  of 
stability  (3.2.8)  and  (3.2.9).  A  continued  investigation  is,  of 
course,  only  called  for  if  these  necessary  conditions  are 
satisfied.  Henceforward  we  may  therefore  assume  in  particular 

P3EU1]  =  0  ,  (3.6.22) 

and  we  may  omit  the  first  term  in  (3.6.20). 

Our  next  step  is  to  show  that  the  quadratic  remainder 
[ V ]  in  (3.6.20)  may  be  made  as  small  as  we  please  in  comparison 
with  P2[v],  by  choosing  the  positive  constants  in  (3.6.3) 
sufficiently  small.  This  is  indeed  the  primary  purpose  of  our 
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reduction  of  the  energy  expression  to  the  form  (3.6.20).  The 
actual  proof  of  our  assertion  is  again  comparatively  simple. 

It  follows  exactly  along  the  same  lines  as  our  proof  in  section 
3.3  that  a  positive  greatest  lower  bound  for  the  quotient  of 
P^Lu]  and  T^Lu]  is  a  sufficient  condition  for  stability.  The 
homogeneous  quadratic  form  satisfies  inequality  (3.3.13)^ 

where  has  been  replaced  by  A^^^  and  y^  by  z^.  The  coefficients 
A\,.j  given  by  (3. 6.1 8)  may  be  made  as  small  in  modulus  as  we 
please  by  choosing  our  positive  constants  in  (3.6.3)  sufficiently 
small.  By  a  suitable  choice  of  these  constants ^  the  absolute 
value  of  X^^z^zM*  may  be  reduced  to  an  arbitrarily  small  positive 
fraction  of  C^^^z^z^j  the  typical  Integrand  of  T^Lv]^  and  we 
oboain  for  the  integrals  in  question  the  inequality 

iRgEvJI  <  s'TgEyj  ,  (3.6.23) 

where  e  *  is  a  positive  constant  as  small  as  v^e  please.  By  means 
of  (3.6.2)  we  obtain  the  desired  result 

|R2[v]i  <  ePglvj  ,  (3.6.24) 

where  e  =  ^'7^2  is  another  positive  constant  as  small  as  we  please. 
Substituting  the  estimate  (3.6.24)  into  (3-6.20),  v;e  have  the 
inequalities 

P[aui+v]  =  a^^P^[u^]  a^P2^[u^,vJ  +  a^P^^[u^,v]  +  a^P2^^[u^,v]  + 

+  (l+E)P2[yJ  +  0(a5)  ,  (3.6.25) 

where  either  both  upper  signs  (  >  and  -  )  or  both  lower  signs 
(<  and  +  )  have  to  be  combined. 
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We  shall  now  investigate  the  functional  of  the  displace¬ 
ment  field  V  orthogonal  to  defined  by 

Pgtvj  +  »  (3.6.26) 

for  a  prescribed  arbitrary  constant  value  of  the  amplitude  a  of 

*) 

the  buckling  mode.  In  particular,  we  shall  determine  the  minimum  ' 
of  (3.8.26)  for  any  given  value  of  a.  Although  the  functionals 
of  V  whj.ch  appear  in  (3.6.25):,  viz. 

(l+e)P2[v]  +  a^P23^[u^,v]  +  v]  , 

(3.6.27) 

are  slightly  different  from  (3.6.26),  we  may  give  a  complete 
description  of  the  behaviour  of  the  functionals  (3.6.27)  in  terms 
of  the  behaviour  of  (3.6.26).  The  sole  difference  between 
(3.8.26)  and  (3.8.27)  is  that  the  latter  functionals  have  an 
additional  factor  (l+e)  in  the  first,  quadratic  term;  the  linear 
terms  are  identical.  If  V  denotes  the  value  of  (3.6.26)  for 
some  displacement  field  v,  then  the  displacement  fields  v/(lTe) 
yield  values  V/(lTe)  for  the  functionals  (3.5.27)^  where  again 
either  the  upper  sign  or  the  lower  sign  has  to  be  applied  every¬ 
where.  Without  any  loss  in  generality  we  may  assume  e  <  1,  since 
e  is  an  arbitrarily  small  positive  number.  The  minima  of  the 
functionals  (3.6.27)  are  then  also  given  by  the  minimum  of  (3.6.26) 
divided  by  (l+s)  . 

in - 

^  Here  again  we  shall  not  prove  the  existence  of  such  a 
minimum,  although  it  will  be  shown  that  a  stationary  value 
of  (3.6.26)  is  necessarily  a  minimum. 
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Let  v'  denote  the  displacement  field  orthogonal  to  the 
buckling  mode  for  which  (3.6.26)  attains  its  minimum.  Let 
denote  an  arbitrary  kinematically  admissible  displacement  field 
subject  to  the  same  condition  of  orthogonality  with  respect  to 
the  buckling  mode 

=  0  •  (3.6.28) 

The  minimum  property  of  (3.6.26)  for  the  displacement  field  v'  Is 
then  expressed  by  the  inequality 

e  +  a2p3iCai:>aJ  +  a^P]ii  Cai.>aJ]  + 

+  e"2p2[fl.]  1  0  ,  (3.6.29) 

holding  for  arbitrary  values  of  the  constant  £ " .  Since  the 
second  term  is  non-negative  for  all  which  satisfy  (3.6.28),  the 
necessary  and  sufficient  condition  resulting  from  this  inequality 
is  that  v=v‘  satisfies  the  equation 

Piita^a)  +  a^P2itHi:-a]  +  a^P3i(Hi'a]  +  a^p^iCai^-a]  =  o 

(3.6.30) 

for  every  kinematically  admissible  displacement  field  which 
satisfies  the  orthogonality  relation  (3.6.28). 

The  existence  of  a  solution  of  the  linear  variational 
equcition  (3.6.30)  will  be  assumed.  It  follows  from  (3.6.29)  that 
this  solution  defines  an  absolute  minimum  of  the  functional 
(3.6.26)  and  is  therefore  necessarily  unique.  Uniqueness  of  the 
solution  of  (3.6.30)  may  of  course  be  confirmed  Indirectly.  Thus 
assume  that  this  equation  has  two  distinct  solutions,  say  v*  and 


60 


v" .  Their  difference  would  satisfy  the  equation 

=  0  .  (3.6.31) 

Taking  t]_  =  (y*-v")  in  (3.6.31)^  we  obtain  by  means  of  (3.^.7) 

2?^[{v^-v") ]  =  0  ,  (3.6.32) 

in  contradiction  to  (3.6.2)  unless  (y*-y")  =  0. 

The  actual  evaluation  of  the  minimum  value  of  (3.6.26), 
once  the  variational  equation  (3.6.30)  has  been  solved,  may  be 
simplified  by  appropriate  use  of  this  equation.  If  we  take  r[ 
equal  to  the  minimizing  displacement  field  y*,  we  obtain  by  means 
of  (3.4.7) 

(3.6.33) 

The  minimum  value  of  (3.6.26)  may  therefore  be  evaluated  by 
either  of  the  following  formulae,  whichever  may  be  most  convenient 
in  a  particular  case 

Min.CPgLv]  H  +  a3p23^[u^,v]  +  = 

=  ■!  {  a'^Pgj^Eu^,  v' ]  +  a^P23^[u^,v' ]  +  a^Pj^j^[u^,  v' ]  ]  = 

=  -Pgtv'J  •  (s.e.s'^) 

The  second  formula  is  also  of  some  general  interest.  It  shows 
that  the  minimum  of  (3.6.26)  is  always  negative  (or  at  least  non¬ 
positive,  if  the  remote  possibility  that  y*  vanishes  identically 
is  also  allowed  for) . 

The  linearity  of  equation  (3.6.30)  implies  that  we  may 
always  write  its  solution  in  the  form 
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v'  =  a^Vg  +  (3.6.35) 

where  the  (vector)  function  is  the  unique  solution  of  the 
equation 

^21 ^  (3.6.36) 

which  satisfies  the  orthogonality  condition  =  0.  We 

omit  the  similar  equations  for  and  Vj^  because  we  shall  not 
need  their  explicit  form.  Substituting  (3.6.35)  into  (3.6.3^) 
we  may  write  the  result  of  our  analysis  of  the  minimum  of  the 
functional  (3.6.26)  in  the  form 

Min.  {  Pglv]  +  +  a^P^j^Euj^^v]  +  a^Pi^^[u^,v]]  = 

=  -  aVg[Vg]  +  0{aP)  ,  (3.6.37) 

where  is  the  unique  solution  of  (3.6.36) . 

We  refrain  from  giving  a  more  explicit  result  for  the 
minimum  (3.6.3?) ^  v/hich  would  be  available  if  full  use  had  been 
made  of  (3.6.34)  and  (3 .6 .35) ^because  we  have  already  disregarded 
in  (3.6.25)  the  explicit  form  of  terms  of  the  fifth  and  higher 
orders  in  the  amplitude  a  of  the  buckling  mode.  Hence  no  useful 
purpose  would  be  served  by  a  more  accurate  evaluation  of  the 
minimum  (3.6.3?).  It  appears  in  fact  that  the  terms  a^P^^[u^^y] 
and  a^P2|^[u^^v]  could  already  have  been  omitted  from  (3.6.25)  and 
(3.6.26)  v/ithout  any  detrimental  effect  on  our  result  (3.6.3?)’. 
The  question  therefore  arises  whether  our  analysis,  leading  from 
(3.6.26)  to  (3.6.3?),  is  perhaps  unnecessarily  complicated  by  the 
retention  of  the  tv/o  terms  under  discussion.  We  did  in  fact 
conjecture  already  at  the  start  of 
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our  analysis  that  only  ’’small"  deviations  from  a  multiple  of  the 
critical  buckling  mode  are  likely  to  be  significant  in  the  dis¬ 
cussion  of  stability..  It  was  for  this  reason  that  an  arbitrary 
kinematically  admissible  displacement  field  u  was  immediately 
written  in  the  form  (3.6.1).  Our  surmise  that  we  may  most  likely 
confine  our  attention  to  "small"  deviations  v  from  a  multiple  of 
the  buckling  mode,  was  also  supported  by  the  appearance  in 
(3.6.25),  in  addition  to  a  quadratic  term  in  y,  of  only  linear 
terms  in  y  multiplied  by  a  quadratic  factor  in  a  or  a  higher  order 

factor  in  a.  The  most  critical  part  of  y  was  therefore  likely  to 
2 

be  of  order  a  ,  a  conjecture  that  is  confirmed  by  our  analysis. 

Nevertheless  we  emphasize  that  all  these  conjectures  are  in  need 

of  adequate  proof .  Such  a  proof  can  only  be  obtained  from  the 

analysis  which  v;e  have  given,  or  from  a  similar  discussion.  It 

3  r  1 

does  not  seem  possible  to  rejeo  .e  terms  a  P^^^Lu^^yJ  and 
a^Pj^l[Ui,y]  in  (3.6.25)  straight  away  without  exposing  the  entire 
argument  to  legitimate  objections . 

V/e  return  now  to  a  final  discussion  of  (3.6.25). 
Remembering  that  the  minimum  values  of  (3.6.27)  may  be  obtained 
from  the  result  (3.6.37)  for  the  minimum  of  (3.6.26)  by  dividing 
by  (l+e),  we  have  the  inequalities 

P[au^+v]  2  ^  (3.6.38) 

Min.  P[au,+v]  <  a^[Pj,[u,  ]  -  Pp[vp]}  +  0(a5)  ,  (3.6.39) 
(a=const.)  ^  ^  ^ 

where  e  is  a  positive  number  that  may  be  made  as  small  as  we 

please  by  selecting  the  positive  constants  A, g,g  *  ,g’’,  . .  .  in 

(3.6.3)  sufficiently  small.  Introducing  the  constant  (cf  (3.6.3^)) 
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=  ^4 fill]  ■**  '*'  ^2^-2^  "  (3.6.40) 

we  obtain  from  (3.6.39)  the  necessary  condition  for  stability  of 
the  critical  state  of  neutral  equlllbrluni 

2  0  ^  (3.6.41) 

and  from  (3.6.38)  the  sufficient  condition  for  stability 

>  0  .  (3.6.42) 

These  statements  are  justified  by  the  fact  that  v;e  may  always 
make  e  so  small  that  the  coefficients  of  a^  in  (3.6.38)  and 
(3.6.39)  have  the  same  sign  as  A^^ .  We  may  then  reduce  the  bound 
A  on  the  modulus  of  a  still  further,  if  necessary,  in  order  to 
ensure  the  domination  of  the  first  term  in  (3.6.38),  (3.6.39) 
over  the  second  term  of  order  a*-^.  It  will  be  clear  that  no 
decision  is  obtained  as  yet  in  the  ambiguous  case  ,  A 

continued  investigation  is  required,  if  it  is  desired  to  obtain 
a  decision  also  in  this  case. 

It  will  be  noticed  that  the  necessary  condition  (3.6.4l) 
requires  more  than  a  positive  definite  fourth  variation  P^i^Lut]. 

It  will  be  remembered  from  the  discussion  in  section  3.2  that  in 
the  case  ^4^^!^  >  0  no  further  necessary  conditions 

could  be  derived  from  (3.2.4).  Condition  (3.6.41)  therefore 
o  contains  the  proof  of  our  previous  statement  in  section  3.2  that 

equilibrium  may  very  well  be  unstable  in  spite  of  the  satisfaction 
of  (3.2.4)  for  every  kinematically  admissible  displacement  field. 
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The  foregoing  derivation  and  proof  of  the  stability 
conditions  (3.6*4l)  and  (3.6.42)  is  based  on  an  auxiliary 
quadratic  functional  T^Eu]  whose  typical  integrand  is  a  positive 
definite  quadratic  form  in  all  arguments  appearing  in  the  typical 
integrand  of  the  energy  Increase  functional  P[u]  defined  by 
(3.1.1).  The  discussion  in  section  3.4,  however,  enables  us  to 
relax  the  restrictions  imposed  on  the  auxiliary  fimctlonal .  Our 
entire  argument  remains  valid,  if  we  replace  T^Eu]  by  any 

*  r  T 

positive  definite  quadratic  functional  T^luJ,  whose  integrands 
are  not  necessarily  positive  definite,  provided  that  we  still 
assume  the  existence  of  a  solution  to  the  original  minimum 
problems  (3.5.1)  and  (3.5.6),  (3.3.7).  The  solution 


* 


u 


1 


Min.  -| - 


0  , 


(3.6.43) 


attained  for  the  displacement  field  u^,  then  also  ensures  a  zero 

solution  of  (3.5.1)  for  the  same  displacement  field  u^  =  u.^^ . 

* 

Moreover,  a  positive  solution  ^  minimum  problem 


* 

=  Min. 


under  the  condition 

T*^[u-^>u]  =  0  , 


(3.6.44) 


(3.6.45) 


also  ensures  a  positive  solution  >  0  of  the  minimum  problem 
(3.5.6),  (3.5.7)^  and  vice  versa. 
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Replacing  (3. 6.1)  by 

u  =  a  +  V  ,  where  3  =  0  (3.6.46) 

all  steps  In  the  previous  analysis  may  be  retraced  up  to  the 
discussion  of  the  remainder  (3.6.19).  Here  we  obtain  with  our 
original  functional  T^Lu]  the  estimate 

.  (3.6.47) 

similar  to  (3. 6.23) ^  and  from  our  original  inequality  (3.6.2) 

Ir^Ev*]!  ie*P2[v*]  ,  (3.6.48) 

‘X’  -M-  f 

where  e  =e  is  a  positive  number  as  small  as  we  please.  Prom 

here  we  may  again  retrace  the  steps  of  our  previous  analysis  . 

* 

Our  final  result  is  now  obtained  in  terms  of  a  constant  , 
defined  by 

H  =  '  (3.6.49) 

•K- 

where  v^  is  the  unique  solution  of  the  equation 

°  (3.6.50) 

which  satisfies  the  orthogonality  condition  T^j[u^,»y  ]  =  0. 

The  restriction  (3.6.28)  on  the  otherwise  arbitrary  kinematically 
admissible  displacement  field  r]_  is  here  of  course  replaced  by  the 
analogous  condition  for  ^ 

=  0  .  (3.6.51) 

The  necessary  condition  for  stability  is  now  Aji^  >_  0,  a  sufficient 
* 

condition  is  Aj^  >  0. 
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It  goes  without  saying  that  the  verdict  on  stability 
must  be  independent  of  the  particular  choice  of  auxiliary 
positive  definite  functional  T^Lu].  Hence  we  must  have  the 
identity 

aJ  =  .  (3.6.52) 

An  analytical  proof  of  this  identity  is  not  immediately 

obvious.  The  different  restrictions  (3.6.28)  and  (3.6.51)^ 

* 

imposed  on  and  ^  j  prevent  an  immediate  comparison  of  the 

solution  Vg  of  (3.6.36)  and  the  solution  Vg  of  (3.6.50).  Removal 

of  these  restrictions  is  also  desirable  for  other  reasons .  We 

shall  prove  therefore  that  these  equations  remain  valid ^  if 
* 

and  are  replaced  by  an  arbitrary  kinematically  admissible 
displacement  field  C  which  is  not  subject  to  any  orthogonality 
condition  of  type  (3.6.28)  or  (3.6.51).  It  will  be  sufficient 
to  confine  our  attention  to  (3.6.36)  because  the  argument  in  the 
case  of  equation  (3.6.50)  is  entirely  similar.  Writing  the 
arbitrary  kinematically  admissible  displacement  field  C  in  the 
form  (cf .(3.^.34)) 

i  =  -I  +  il  =  (3.6.53) 

where  j]  satisfies  (3.6.28),  we  may  evaluate  the  left-hand  member 
of  (3.6.36)  for  its  solution  v^,  and  after  replacement  of  j]  by  C 

=  P21^-1"-3^  •  (3.6.5^) 

The  first  term  within  the  braces  vanishes  because  u,  is  the  buck- 

— ± 

ling  mode  (cf .  (3 .5 .3) ) :»  and  the  second  term  within  these  braces 


is  also  zero  on  account  of  (3-6.22) 

"  3P3[ui]  =  0  (3.6.55) 

'Che  last  two  terms  cancel  because  Vg  is  the  solution  of  (3-6.36). 
Hence  we  have  proved  that  Is  also  a  solution  of  the  modified 
variational  equation 

+  P23_[u3_,C]  =  0  ,  (3.6.56) 

holding  for  every  kinematically  admissible  displacement  field  C. 

In  fact,  v^  is  the  unique  solution  of  (3-6.56)  which  also  satisfies 
the  orthogonality  relation 

TiiCui^Vg]  =  0  .  (3.6.57) 

Likewise,  v^  is  the  unique  solution  of  the  same  equation  (3-6.56) 
which  also  satisfies  the  orthogonality  relation 

T*i[ui,V2j  =  0  .  (3-6.58) 

Equation  (3-5-56),  satisfied  by  both  v^  and  v^,  enables 
us  to  compare  these  solutions .  Their  difference  must  be  a  solu¬ 
tion  of  the  equation 

=  0  .  (3.6.59) 

Since  the  critical  buckling  mode  is  simple,  the  general  solution 

* 

for  the  difference  of  v^  and  v^  is 

Vg-Vg  =  bu^  ,  (3.6.60) 


where  b  is  an  arbitrary  constant .  Hence  v^e  have 


if  again  appropriate  use  is  made  of  (3.5*3)  and  of  (3.5.2). 
Substitution  of  ( 3.6.61)  into  (3.6.49)  completes  our  proof  of 
(3.6.52) . 

A  second  reason  why  the  modified  variational  equation 
(3.6.56)  is  to  be  preferred  to  the  original  equation  (3.6.36)  or 
(3.6.50)  is  that  the  modified  equation  is  equivalent  to  a  set  of 
differential  equations  and  boundary  conditions .  These  differential 
equations  and  boundary  conditions  are  nearly  always  indispensable 
in  the  solution  of  any  particular  problem  of  stability  in  a 
critical  case  of  neutral  equilibrium.  Their  derivation  from 
(3.6.56)  proceeds  along  the  same  lines  as  the  derivation  of  the 
equations  of  neutral  equilibrium  from  the  variational  equation 
(3.5.3).  The  fact  that  the  first  term  in  (3.6.56)  and  the  left- 
hand  member  in  (3.5.3)  are  identical  implies  that  the  unknown 
displacement  field  v  and  its  derivatives  appear  in  our  present 
equations  in  exactly  the  same  form  as  u  and  its  derivatives  appear 
in  the  equations  of  neutral  equilibri\im.  The  only  difference  is 
that  the  present  equations  are  no  longer  homogeneous . 

A  curious  difficulty  is  apparently  raised  here  by  the 
general  theory  of  linear  differential  equations.  The  homogeneous 
equations  have  a  non-vanishing  solution  in  the  form  of  the 
critical  buckling  mode  because  they  are  identical  with  the  equa¬ 
tions  of  neutral  equilibrium.  In  order  that  a  solution  of  the 
non-homogeneous  equations  may  exist,  their  prescribed  right-hand 
members  have  then  to  satisfy  a  certain  necessary  condition  in  the 
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form  of  an  Integral  relation.  '  The  question  then  arises  whether 
this  relation  is  actually  satisfied  by  the  right-hand  members  of 
the  equations  resulting  from  the  variational  equation  (3.6.56), 

The  answer  to  the  latter  question  is  fortunately 
affirmative j  and  the  difficulty  raised  by  the  general  theor^r  of 
linear  differential  equations  is  therefore  no  more  than  apparent . 
In  the  interest  of  future  applications,  it  will  be  convenient  to 
prove  this  statement  in  a  slightly  more  general  context  than  is 
actually  needed  for  the  purpose  at  hand.  Our  variational  equation 
(3.6.56)  is  a  particular  case  of  the  more  general  equation 

=  0  ,  (3.6.62) 

where  Is  a  linear  functional  of  the  arbitrary  kinematically 

admissible  displacement  field  X,*  associated  non-homogeneous 

linear  differential  equations  and  boundary  conditions  admit  a 
solution,  only  if  the  right-hand  members  satisfy  a  single  Integral 
relation.  Since  the  right-hand  members  result  from  the  term 

in  (3.6.62),  it  follov73  that  a  solution  exists,  only  if  the 
functional  M^[C]  obeys  an  equivalent  condition.  The  explicit  form 
of  the  latter  condition,  however,  is  obtained  Immediately  from 
(3.6.62)  by  observing  that  this  equation  must  also  hold  for  ^=u^ . 
The  first  term  in  (3.6.62)  vanishes  in  this  case,  and  it  follows 

in - 

^This  necessary  condition  is  usually  called  an  orthogonality 
condition  in  the  theory  of  differential  equations.  We  shall 
avoid  this  terminology  in  order  to  preventthe  danger  of  con¬ 
fusion  with  our  concept  of  orthogonality  based  on  the  auxiliary 
quadratic  functional  T^Luj. 
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that  the  linear  functional  must  satisfy  the  condition 

0  .  (3.6.63) 

in  order  that  the  variational  equation  (3.6.62)  may  have  a  solu¬ 
tion.  Condition  (3.6.63)  must  therefore  be  the  equivalent  of  the 
integral  relation  to  be  satisfied  by  the  right-hand  members  of 
the  differential  equation  and  bouindary  conditions .  It  is  easily 
verified  by  means  of  (3.6.55)  that  in  the  particular  case  under 
consideration, 

,  (3.6.64) 

condition  (3.6.63)  is  indeed  satisfied.  It  follov;s  that  the 
right-hand  members  of  the  differential  equations  and  boundary 
conditions  associated  v/ith  (3.6.56)  actually  satisfy  the  necessary 
condition  in  question  for  the  existence  of  a  solution  of  the  non- 
homogeneous  equations . 

There  is  a  deeper  reason,  apart  from  the  vanishing  of 
in  the  particular  case  under  consideration,  why  the  non- 
homogeneous  variational  equations  of  the  type  (3.6.62),  which  we 
shall  encounter  in  various  special  forms  in  the  sequel,  always 
satisfy  the  Integral  relation  which  is  a  necessary  condition  for 
the  existence  of  a  solution  of  the  associated  non-homogeneous 
differential  equations  and  boundary  conditions  .  It  may  therefore 
be  worthwhile  to  pursue  the  discussion  a  little  further,  even  if 
this  is  not  strictly  necessary  for  our  immediate  purposes . 

Variational  equations  of  type  (3.6.62)  always  enter  our 
discussion  as  the  necessary  and  sufficient  condition  for  the 
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existence  of  a  solution  of  a  minimum  problem 

PgCv]  +  M^Lv]  =  Min.  (3.6.65) 

under  the  side  condition 

Tli[Ui.v]  =  0  ,  (3.6.66) 

where  lVI^[v]  is  some  linear  functional  of  v.  We  shall  assume  that 
its  typical  integrand  depends  only  (and  that  of  course  in  linear 
form)  on  the  displacement  components  and  their  derivatives  which 
also  occur  in  the  energy  increase  functional  (3.I.I)  and  in  T^Eu]. 
It  is  easily  verified  that  the  minimum  problem  (3.6.65),  (3.6.66) 
is  properly  formulated  in  the  sense  that  the  existence  of  a  lower 
bound  of  (3.6.65)  is  ensured.  In  fact,  we  have 

?^lv]  +  M^[v]  >  +  M^[v]  ,  (3.6.67) 

and  the  typical  integrand  of  the  right-hand  member  in  (3.6.67)  is 
obviously  bounded  below.  Let  v'  denote  the  solution  of  (3.6.65)^ 
whose  existence  will  now  again  be  assunied,  and  let  again  denote 
an  arbitrary  kinematically  admissible  displacement  field  which 
satisfies  the  orthogonality  condition  (3.6.28).  The  minimum 
property  of  v’  is  then  expressed  by  the  inequality  of  well-known 
type 

e{Pii[v'Jl]  +  >  0  ,  (3.6.68) 

where  e  is  an  arbitrary  constant.  The  necessary  and  sufficient 
condition  resulting  from  this  inequality  is  the  equation 


(3.6.69) 
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satisfied  by  v=v*  for  every  displacement  field  ji  subject  to 
restriction  (3.6.28).  In  order  to  obtain  from  (3.6.69)  a  varia¬ 
tional  equation  In  which  the  arbitrary  kinematically  admissible 
displacement  field  ^  Is  not  subject  to  the  restriction  (3.6.28) ^ 
we  substitute  (3.6.53)  in  the  left-hand  member  of  (3.6.69),  and 
evaluate  the  result .  We  obtain 

+  +  M^[2i]  =  -I  , 

(3.6.70) 

where  we  have  used  (3.5.3)  and  (3.6.69).  Writing  this  result  in 
the  form 

-  I  =  0  ,  (3.6.71) 

we  have  obtained  an  equation  of  type  (3.6.62),  in  which 

=  M^[^]  -  i  (3.6,72) 

automatically  satisfies  condition  (3.6.63).  We  need  not  fear 
therefore  that  the  right-hand  members  of  the  differential  equa¬ 
tions  and  boundary  conditions  associated  with  (3.6.71)  would 
violate  the  integral  relation  which  is  a  necessary  condition  for 
the  existence  of  a  solution.  The  satisfaction  of  this  condition 
is  in  fact  a  consequence  of  the  proper  formulation  of  minimum 
problem  (3.6.65).  (3.6.66). 

Finally,  a  word  of  caution  may  be  in  order  here  to 
prevent  the  misunderstanding  that  the  foregoing  discussion  proves 
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the  existence  of  a  solution  of  minimum  problems  of  type  {3.6.65) j 
(3.6.66),  their  associated  variational  equations  of  type  (3.6.62) 
or  the  corresponding  differential  equations  and  boundary  condi¬ 
tions.  The  quite  difficult  issue  of  existence  of  solutions  is 
far  beyond  the  scope  of  our  discussion.  Our  purpose  has  merely 
been  to  show  that  the  well-known  necessary  conditions  for  the 
existence  of  a  solution  of  non -homogeneous  differential  equations 
and  boundary  conditions,  are  always  satisfied  in  our  equations 
resulting  from  properly  formulated  minimum  problems . 

The  actual  application  of  the  stability  criterion, 
specified  by  the  constant  (3.6.4o),  is  extrejmely  simple  in  the 
case  of  our  first  example,  the  bar  under  compressive  end  loads 
(fig.  3.1)-  A  critical  case  of  neutral  equilibrium  occurs,  if 
and  only  if  the  compressive  load  N  equals  the  Euler  load 

N  =  N  =  ’4^  .  (3.6.73) 

The  corresponding  buckling  mode  (3.5.20)  for  k=l  may  be  normalized 
here  by  taking  the  maximum  deflection  equal  to  unity 

w^  =  sin  ux/^  .  (3.6.74) 

The  third  variation  vanishes  Identically  in  this  example 

(cf.  (3.2.12)).  Hence  the  second  term  in  the  variational  equation 
(3.6.56)  is  also  identically  zero.  The  general  solution  of  this 
equation  is  a  multiple  of  the  buckling  mode  (3-6.74),  and  it 
follov;s  that  v^  =  0  in  view  of  the  required  orthogonality  to  the 
buckling  mode .  flence  the  constant  is  in  this  case  supplied 
immediately  by  the  fourth  variation  P|^[w],  evaluated  for  the 
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buckling  mode  (3.6.74).  We  obtain 


^4  =  ^ 


tO  P 

B  ^.2  Tix  ^^2  Tix 
L — ^ —  sin  —j^  cos 

O 


T 


^  cos  ^^Jdx  — 


(3.6.75) 


64.0^ 

This  constant  A^^  is  positive,  and  the  critical  case  of  neutral 
equilibrium  at  the  Euler  load  is  therefore  stable. 

The  second  example,  the  three-dimensional  elastic  body, 
is  far  more  complicated.  No  general  solution  of  the  equations  of 
neutral  equllibriiom  is  available  in  this  case .  If  the  buckling 
mode  in  a  critical  case  of  neutral  equilibrium  is  again  denoted 
by  u^,  the  displacement  components  by  u^^  (i=l,2,3)^  it  would  of 
course  be  possible  to  derive  the  differential  equations  corre¬ 
sponding  to  the  variational  equation  (3.6.56)  along  the  same  lines 
as  the  derivation  of  the  equations  of  neutral  equilibrium  in 
section  3.5.  fbe  resulting  equations  are  extremely  complicated, 
and  little  additional  insight  is  gained  from  them.  Let  it  there¬ 
fore  suffice  to  give  the  appropriate  expression  for  the  second 
term  in  (3.6.56).  The  interested  reader  may  then  derive  the 
associated  explicit  differential  equations  and  boundary  conditions 
From  (3.2.18)  we  obtain  in  this  example 


+ 


l-2v  ^lk,k'^lh,i*=h,  j  l-2v  ^Ih,  • 


(3.6.76) 
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A  final  remark  seems  appropriate  before  closing  the 
present  section.  It  should  be  noted  that  the  detailed  investiga¬ 
tion  of  stability  in  a  critical  case  of  neutral  equilibrium  is 
necessary  only  in  cases  which  ought  perhaps  to  be  classified  as 
"exceptional"  from  a  mathematical  point  of  view.  The  "general" 
case  in  a  mathematical  sense  would,  of  course,  be  that  the  third 
variation  ^oes  not  vanish  for  the  buckling  mode  u^,  that  is 

P3[u^]  /  0  .  (3.6.77) 

Instability  of  the  critical  state  of  neutral  equilibrium  is  then 
immediately  obvious  from  the  necessary  condition  for  stability 
(3.2.7)  or  (3.6.22),  and  any  further  investigation  is  superfluous. 
Nevertheless,  many  problems  do  exist  in  which  this  necessary 
condition  (3.2.7)  happens  to  be  satisfied,  and  the  continued 
analysis  of  the  present  section  is  indeed  required  in  such  cases . 
The  vanishing  of  is  usually  caused  by  certain  symmetry 

properties  of  the  structure  and/or  its  buckling  mode,  which  over¬ 
rule  the  prediction  that  (3.6.77)  is  the  general  case  in  a 


mathematical  sense. 
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3 . 7  The  Indecisive  case. 

In  the  preceding  section  we  have  derived  necessary  and 

sufficient  conditions  for  the  stability  of  the  fundamental  state 

in  a  critical  case  of  neutral  equilibrium  characterized  by  a 

simple  critical  buckling  mode  u^ .  Tlie  discussion  of  the  case  of 

multiple  buckling  modes  is  postponed  until  the  next  section.  The 

necessary  and  sufficient  conditions  were  obtained  in  the  form  of 

the  Inequalities  ^  0  and  >  0  respectively^  where  is  a 

constant  defined  by  (3.6.4o).  We  shall  now  examine  how  the 

investigation  may  be  continued  in  the  indecisive  case  A^=0. 

At  a  first  glance  it  may  appear  sufficient  for  this 

purpose  to  refine  the  accuracy  of  the  estimates  (3.6.9)  and 

(3.6.37)  by  taking  into  account  explicitly  the  terms  of  order 
5 

a"^ .  It  would  then  be  overlooked^  however,  that  the  remainder 
(3.6.19)  in  (3.6.20)  also  involves  a  term  of  order  a^^  if  the 
solution  (3.6.35)  is  substituted  into  (3.6.I9).  A  more  careful 

o 

analysis  is  therefore  required.  We  may  of  course  expect  that  the 
decision  on  stability  will  be  given  by  a  displacement  field  close 
to 

u  =  au^+a^Vg  ,  (3.7.1) 

which  barely  fails  to  yield  the  required  verdict .  Hence  it  will 
obviously  be  expedient  to  write  an  arbitrary  kinematically 
admissible  displacement  field  in  the  form 

u  =  au^+a  v^+w  ,  v/here  T^^[u^,w]  =  0  (3.7.2) 
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We  need  a  convenient  notation  for  the  expansion  of  a 
functional  S^[u]j  whose  integrand  is  a  homogeneous  polynomial  of 
degree  m  in  the  displacement  components  u^  and  their  derivatives^ 
similar  to  (3.^.6)^  but  now  for  the  sum  of  three  displacement 
fields.  We  write 


S  [u+y+w]  =  Z 

(p+q+r=m) 


pqr 


[u, v,w]  , 


(3.7.3) 


where  the  typical  Integrand  of  Sp^^Eu.Vjw]  contains  all  terms  of 
degree  p  in  u^  of  degree  q  in  y>  and  of  degree  r  in  w.  Vie  note 
the  six  symmetry  relations 


S 

S 


pqr 


prq 


[u,w,v] 


=  Sqp^[v,u,w]  =  S 
=  Sq^p[v,w,u]  =  S 


rqp 


[w,v,u]  , 


and  the  Identities  similar  to  (3.^.7) 


(3.7.^) 


Soqr[u,v,w] 

Spqr[:ii^H>w] 


(3.7.5) 


p^q?^  ’  ^(p+q)r^^^-^'  “  p!q!r! 


(3.7.6) 


The  arguments  of  the  typical  Integrand  F(2^)  of 
P[u]  are  written  in  a  form  corresponding  to  (3.7.2)  and  similar 
to  (3.6.6) 


=  ay 


^a^z^-s^  ,  or  ^  =  a;i^j^+a^Zg+s  ,  (3.7.7) 


vrhere  now  represents  the  arguments  corresponding  to  y^ .  The 
Taylor-expansion  of  the  typical  Integrand  is  now  written  in  the 


form 
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+  i  s^sH-P^^^ fa^j^+a^g)  'g  ’  (3.7.8) 

where  the  arguments  in  the  third  order  derivative  are  given  by 

ay^+a^Z2+¥^  =  ay^+a^Zg+G^'s^  ,  (X  not  summed;  0  <  0^  ^  l)  , 

(3.7.9) 

The  second  derivatives  in  (3.7.8)  are  rev/rltten  as 

F^^^(aii+a2z2)  =  P^^^(O)  +  (ay^+a^zpP^^^^  (Oa^i+Oa^)  , 

(0  <  0  <  1)  .  (3.7.10) 

V/e  may  now  proceed  to  identify  the  various  terms  after 
the  collection  of  all  integrals .  V/e  obtain  the  symbolic  relations 

^(aij^+a^Zg)  =:^  a^P^Euj^]  +  a^{P|^[u^]  +  P23^[u^»V2]  +  ?p[Vp]}  + 

+  a5{P^[u^]  +  +  Pj^gEu^^Vg]]  + 

+  a^{Pg[u^]  +  Pjj^^Eu^jVg]  +  P22[u2^l2^  ^3(^2]}  + 

+  0{a7)  .  (3.7.11) 

s^P^j^(a£^+a^Z2)  =i>  a^EPg^^Eu^,^  ]  +  P^^j^Evg^w]]  + 

+  a^EP^jEu^.w]  +  P3^]^]^Eu2>l2'-b  + 

+  a^{P2^j^Eu^,w]  +  PgiiEuj^.Vg.w]  +  Pgi^EVg,^]}  + 

+  a^{P^j^Eu^,w]  +  Pjj^j^Eu^, Vg'-^  Pl21^-l'-2"-b  ■*■ 

P4ll[2il'Z2'l'3  +  P221^-1'-2"-^  ^31(22^^^  + 

+  0(a'^)  ,  (3.7.12) 
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i  =^P2[w]  , 

For  the  remainder  of  (3.7-8)  we  have  the  relation 


.V  Rglw]  , 


where  Is  defined  by 

+  5  (ayi+a^z^)P^^^^(©a;ir^+ea%)  . 


(3.7.13) 


(3.7.1^^) 


(3.7.15) 


Before  we  write  down  our  resulting  expression  for  the 
energy  increase  in  the  displacement  field  (3.7.2)  from  the  funda¬ 
mental  state  to  some  state  II,  it  is  convenient  to  make  use  of 
the  simplifications  available  from  our  previous  investigation. 

We  observe  that  the  first  term  and  the  second  term  in  (3.7.11) 
vanish  by  virtue  of  (3.6.22)  and  the  fact  that  in  our  indecisive 
case 

^21  ^-1^-2 ^  =  0  .  (3.7.16) 

Restricting  our  attention  first  to  a  displacement  field  (3*7.1) ^ 
i.e.  to  w  =  0,  the  entire  energy  increase  is  given  by  (3.7.11). 

A  necessary  condition  for  stability  is  then  that  the  lowest 
power  of  a  in  (3.7.11)  is  even  and  has  a  non-negative  coefficient. 
Hence  we  obtain  the  necessary  conditions  for  stability  in  the 
indecisive  case  A)j=Q 

P5[HiJ  +  °  ’  (3.7.17) 

+  P22^-1"-2^  ^3l--2^  =  °  •  (3.7.18) 
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Henceforward  we  shall  assume  that  these  necessary  conditions  are 
satisfied. 

If  we  now  consider  (3«7.12)j  we  observe  that  the  first 
term  vanishes  identically  because  is  the  solution  of  (3.6.36). 
Moreover. »  our  previous  experience  with  the  terms  v/hich  are  linear 
in  v  in  (3.6.20)  indicates  that  only  the  term  with  the  lowest 
power  of  a  as  a  factor  is  actually  significant.  We  shall  profit 
from  this  experience  by  omitting  in  (3.7.12)  all  terms  except  the 
first  non-vanishing  term.  We  emphasize^  however,  that  this  omis¬ 
sion  lacks  a  proper  Justification  at  this  stage .  A  rigorous 
analysis  would  require  retention  of  all  terms  in  (3.7.12),  and 
it  would  then  appear  a  posteriori  that  only  the  term  with  a*^  as 
a  factor  is  significant.  It  may  be  left  to  the  reader  to  complete 
the  proof  in  this  sense. 

Our  resulting  simplified  expression  for  the  energy 
increase  (3.1.1)  in  the  displacement  field  (3.7.2)  reads 

P[au^-i-a^V2+w]  =  a^{Pg[u^]  +  Pi^^lu^^v^]  + 

+  a^{P23^[u^,w]  + 

+  Pgtw]  +  •  (3.7.19) 

The  further  analysis  is  the  complete  analogue  of  the  Investigation 
of  (3.6.20)  in  the  preceding  section,  and  it  may  again  be  left  to 
the  reader  to  fill  in  the  details  of  the  following  discussion. 

It  is  first  established  that  the  remainder  in  (3.7.19)  satisfies 
the  inequality 


81 


|R2[w]|  1  ePgtw]  ,  (3.7.20) 

where  e  is  a  positive  mimber  as  small  as  v;e  please. 

Let  denote  the  solution  of  the  minimum  problem 

PgEw]  +  ^  Min.,  (3.7.21) 

where  the  displacement  field  w  is  restricted  by  the  side  condition 

T^^[u^,w]  =  0  .  ■  (3.7.22) 

The  displacement  field  is  a  solution  of  the  variational  equa¬ 
tion 

=  0  ,  (3.7-23) 

holding  for  every  kinematically  admissible  displacement  field 
subject  to  the  orthogonality  restriction  (3.6.28).  The  solution 
is  the  unique  solution  of  (3.7.23)  which  also  satisfies  the 
orthogonality  condition  (3.7.22) .  The  variational  equation 
(3.7.23)  may  again  be  replaced  by  a  similar  equation  in  terms  of 
an  arbitrary  kinematically  admissible  displacement  field  ^  which 
is  not  subject  to  the  orthogonality  restriction  (3-6.28)  on  r^. 

In  view  of  the  identity 

^  +  2P2j^Lu^,V2]  = 

=  -  P2t'^2h  =  0  ^  (3.7.24) 

resulting  from  (3-6.40)  and  Aj^=0,  the  latter  equation  has  exactly 
the  same  form  as  (3-7-23) 

+  P31CH1.L)  +  =  0 


(3.7.25) 
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Introducing  the  constant  defined  by 

Ag  =  '''  ^22^-1 '-2^  ^^3^-2^  ■  ^2^-3^  ' 

(3.7.26) 

a  further  necessary  condition  for  stability  in  the  Indecisive 
case  A2|^=0  is 

Agio,  (3.7.27) 

a;ad  a  suff Icleint  condition  for  stability  is 

I 

Ag  >  0  .  (3.7.28) 

No  decision  is  obta.ined  once  again;,  if  the  constant 
A^  equals  zero,  and  the  analysis  has  to  be  continued  anew  in  this 
indecisive  case  of  higher  oi*der.  The  procedure  to  be  followed 
then  will  be  clear  from  the  foregoing  analysis.  It  does  not 
present  any  new  difficulties,  although  the  energy  expression  for 
a  displacement  field 

u  =.  au^+a  v^+a'^w.^+w  ,  where  =  0  ,  (3.7.29) 

and  the  corresponding  equations  become  increasingly  more  complex. 
We  may  omit  a  further  discussion,  the  more  so  since  the  occurrence 
of  this  quite  exceptional  case  seems  highly  unlikely  in  practical 
applications  of  the  theory. 
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3.8  The  case  of  multiple  buckling  modes. 

The  investigation  of  stability  in  a  critical  case  of 
neutral  equilibrium  in  sections  3.6  and  3.7  has  been  restricted 
to  the  case  of  a  simple  critical  buckling  mode.  We  shall  now 
discuss  the  slightly  modified  analysis  which  is  required  in  the 
case  of  multiple  buckling  modes.  Let  n  denote  the  number  of 
linearly  independent  critical  buckling  modes .  The  first  n  solu¬ 
tions  in  the  sequence  of  minimum  problems  (3.^.26)j  (3. ^.2?)  are 
then  all  zero,  and  the  solution  of  the  (nH'l)-st  problem  is 


positive 


tj  =0 

n 


(n+l) 


(3.8.1) 


The  buckling  modes  u^^  (h=l,2,...n)  are  the  minimizing  displace- 
ment  fields  in  the  first  n  problems,  noraalized  by  the  condition 

It  is  expedient  here  to  represent  an  arbitrary  kine¬ 
matically  admissible  displacement  field  u  by  means  of  (3. ^.23) 
in  the  fom 

u  =  a^j^+v  ,  where  T^^[uj^,v]  =  0  ,  h=l,2,...n  (3.8.2) 

Here  and  in  the  remainder  of  this  section  we  shall  employ  the 
summation  convention  for  a  repeated  Latin  subscript  in  the  sense 
that  such  a  repeated  subscript  implies  summation  from  1  to  n^ 
in  the  absence  of  an  explicit  statement  to  the  contrary. 

The  arguments  y^  of  the  typical  integrand  F{^)  of  P[u], 
defined  in  section  3.2,  are  written  in  a  form  corresponding  to 
(3.8.2),  and  similar  to  (3.6.6) 
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=  aj^y^  +  ,  or  +  z  ,  (3.8.3) 

where  represents  these  arguments  in  the  buckling  mode  Uj^. 

The  analysis  now  follows  again  along  the  same  lines  as  in  section 
3.6  with  the  single  exception  that  au^  and  a;^^  are  replaced  every¬ 
where  by  the  sums  and  a^j^.  The  first  important  result  is 

an  expression  for  the  energy  increase  (3.1.1)^  similar  to  (3.6.20) 

P[ajUj^+v]  =  P3[ajU^]  +  P4[aj^u^]  +  +  P3I  v]  + 

+  P4i[aj^j^jv]  +  P2[v]  +  RgEv]  +  O(a^)  ,  (3.8.4) 

v;here  the  symbol  O(a^)  stands  for  a  term  which  tends  to  zero  as 
a  homogeneous  polynomial  of  the  fifth  degree  in  the  amplitudes 
a.j^  of  the  buckling  modes,  if  these  amplitudes  themselves  approach- 
zero.  An  Inequality  similar  to  (3.6.24)  may  be  proved  for  the 
term  R2[v],  viz. 

iRgtvJI  S  EPgfv)  «  (3.8.5) 

where  e  is  again  a  positive  number  as  small  as  we  please.  In 
the  present  case  this  inequality  is  based  on  the  inequality  for 

“(n+l)  >  0 

PgCv]  2  “(n+l)'^2^-^  ’  (3.8.6) 

which  replaces  (3.6.2). 

It  is  convenient  to  derive  immediately  some  necessary 
conditions  for  stability  by  considering  a  displacement  field  u 
which  is  a  linear  combination  of  the  critical  buckling  modes, 
i.e.  by  taking  v  s  0  in  (3.8.4).  We  must  obviously  have  then 
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(3.8.7) 

(3.8.8) 


for  every  choice  of  the  amplitudes  We  Introduce  the  con- 

stants  ^ 


A 


hij 


■5  ^  ® 


hi  jk 


1 


^llll^-h-’-i-’-j^-k^  "  (3.8.9) 


where  the  last  functional  is  defined  as  the  sum  of  all  terms  in 
the  expansion  of  P]^[u]  for  li=u^+Ujl^+Uj+Uj^j  which  are  linear  in  all 
four  constituents  Uj^^u^.?Hj  Condition  (3.8.?)  is  then 

expressed  by  the  requirement  that  the  cubic  form 


A,  .  .a,  a. a ,  =  0 
hlj  h  1  J 


(3.8.10) 


must  vanish  identically .  Hence  we  must  have 

Aj^^j  —  ^111  h^~i J  ^  ^  ...n  »  (3*8,11) 

Likewise j  condition  (3.8.8)  requires  that  the  quartlc  form 

^  °  (3.8.12) 


is  non-negative  for  every  combination  of  amplitudes  of  the  buck¬ 
ling  modes .  Henceforward  v;e  shall  assume  that  both  conditions 
(3.8.11)  and  (3.8.12)  are  satisfied. 

Again  similar  to  (3.6.26)  we  now  investigate  the  func¬ 
tional  of  the  displacement  field  v  orthogonal  to  all  buckling 
modes,  defined  by 


’  We  reserve  the  symbol  later  use. 
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(3.8.13) 

for  prescribed  arbitrary  constant  values  of  the  amplitudes  of 
the  buckling  modes.  Let  v*  be  the  displacement  field  v  which 
minimizes  the  functional  (3.3.13) .  It  is  the  unique  solution, 
orthogonal  to  all  critical  buckling  modes,  of  the  variational 
equation 

(3.8.14) 

holding  for  every  kinematically  admissible  displacement  field  £ 
which  satisfies  the  requirement  of  orthogonality  to  all  buckling 
modes 

Tll[Uh.il]  =  0  .  h=l,2,...n  .  (3.8.15) 

Equation  (3.8.14)  Is  again  linear  in  v.  We  write 

=  I  •  (3.8.16) 

The  solution  of  (3.8.14)  may  then  be  written  in  the  form 

^h^i-hi  ’  (3.8.17) 

where  is  the  unique  solution,  orthogonal  to  all  buckling 
modes,  of  the  equation 

+  §  Pxil^-h’%".^^  ""  °  "  (3.8.18) 

and  0(a'^)  stands  for  a  contribution  that  tends  to  zero  for 
a^  — ■>  0  as  a  homogeneous  polynomial  of  the  third  degree  in  a^. 
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Equation  (3,8,l8)  may  again  be  replaced  by  a  similar  equation  in 
terms  of  a  completely  arbitrary  kinematical  displacement  field 
C  which  is  not  subject  to  any  orthogonality  restriction.  In 
view  of  (3.8.11)  this  equation  is  again  identical  in  fom  to 
(3.8.18) 

=  0  •  (3.8.19) 


The  minimum  value  of  (3.8.13)  may  now  be  written  in 
either  of  the  forms 


'^hijkVi^A  ^  (3.8.20) 


where  the  coefficients  C,  .  are  defined  by 

hijk 


ilr  O  1  * 


(3.8.21) 


'hiJk  2  ll"-hi^~jk' 

From  the  inequality  (3.8.5)  we  nov/  obtain  two 
inequalities  similar  to  (3.6.38)  and  (3.6.39) 

^  ^®hijk;  '  ‘^hl jk^^h^i^j®-k  ''' 

(3.8.22) 

(„  ^  ^®hijk  ■  SiJk^^h^i^J^k 

^  (3.8.23) 


where  e  is  again  a  positive  number  as  small  as  we  please.  Intro¬ 
ducing  the  constants 


\ijk  ^hijk  ’  ^hljk  ^ 


(3.8.24) 


we  obtain  from  (3.8.23)  the  necessary  condition  for  stability  of 
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the  critical  state  of  neutral  equilibrium  with  multiple  buckling 
modes 

^  °  (3.8.25) 

for  all  values  of  the  amplitudes.  Likewise^  vie  obtain  from 
(3.8.22)  as  a  sufficient  condition  for  stability 

^hilk^h^l^-j^k  positive  definite.  (3.8.26) 

No  verdict  on  stability  is  again  obtained^  if  the 
quartic  form  in  (3.8.25)  and  (3.8.26)  is  semi -definite .  A 
further  investigation  would  be  required  in  03.'*der  to  get  a  decision. 
Such  an  investigation  is  more  complicated  than  the  investigation 
of  the  indecisive  case  for  a  simple  buckling  mode  because  the 
semi -definite  quartic  form  presimiably  is  zero  only  for  one  or 
more  special  sets  of  values  for  the  ratios  of  the  amplitudes  of 
the  buckling  modes .  V/e  shall  not  pursue  this  investigation  since 
the  indecisive  case  for  multiple  buckling  modes  is  indeed  a  rare 
exception.  In  the  author’s  experience  it  has  never  been 
encountered  in  actual  stability  problems . 

Finally^  it  ;-7ill  be  remembered  from  our  discussion  at 
the  end  of  section  3.6  ^hat  the  vanishing  of  the  third  variation 
ought  to  be  regarded  as  an  "exceptional”  case  in  a  mathematical 
sense.  Hence  it  is  even  more  "exceptional"  in  the  same  sense 
if  all  n(n+l)(n+2)  coefficients  (3.8.11)  vanish.  The 

symmetry  properties  of  the  structure  and/or  its  buckling  modes 
may  of  course  again  overrule  this  consideration.  It  may  again 


89 

result  from  such  symmetry  that  (3.8.11)  is  indeed  satisfied. 
Nevertheless,  experience  indicates  a  more  frequent  occurrence  of 
instability  from  a  violation  of  (3.8.11)  than  from  (3.6.77)  in 
the  case  of  a  simple  buckling  mode. 
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3 . 9  Approximate  solution  of  stability  problems . 

The  analysis  in  the  present  chapter  is  essentially 
mathematical  in  character^  even  if  we  have  employed  the 
terminology  of  the  theory  of  stability.  No  approximative 
physical  assumptions  of  any  kind  have  been  made  so  far.  In  this 
connection  we  may  ignore  the  two  examples,  which  have  been 
discussed  merely  for  the  purpose  of  some  Illustration  of  the 
argument.  In  view  of  the  absence  of  any  approximation,  the 
analysis  is  entirely  rigorous,  except  for  the  omission  of  proofs 
of  existence  of  solutions  for  the  minimum  problems  encountered 
in  sections  3.^  and  3.6.  Moreover,  the  analysis  is  equally 
applicable  to  other  problems  in  the  calculus  of  variations, 
specified  by  some  functional  P[u]  which  takes  the  place  of  our 
energy  increase  functional  (3.1.1). 

It  will  cause  no  surprise,  however,  that  the  actual 
application  of  the  theory  to  problems  of  elastic  stability 
requires  far-reaching  simplifications .  The  number  of  problems 
which  are  capable  of  a  rigorous  solution  is  small  even  in  the 
classical  linear  theory  of  elasticity.  The  limitations  of  the 
available  analytical  tools  become  even  more  apparent  in  the 
more  complicated,  inherently  nonlinear,  theory  of  elastic 
stability.  Approximations  of  some  type,  and  in  most  cases  even 
of  several  types,  are  nearly  always  Indispensable  in  order  to 
obtain  concrete  answers  to  the  stability  question  in  specific 
problems.  It  seems  therefore  appropriate  to  end  the  present 
chapter  with  a  preliminary  discussion  of  some  of  the  more 
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important  types  of  approximations  employed  in  the  theoi'y  of 
elastic  stability.  This  general  discussion  will  be  helpful  in 
our  future  choice  of  suitable  approximations  in  particular  cases, 
and  it  will  shield  us  to  some  extent  from  incautious  approaches. 

It  cannot  be  avoided  here  that  parts  of  the  argument  are  perhaps 
less  meaningful  to  readers  without  any  previous  experience  in 
the  theory  of  elastic  stability.  It  may  be  some  comfort  to 
such  readers  that  we  shall  repeatedly  return  to  these  questions 
in  future  chapters,  where  the  details  of  the  analysis  will 
clarify  the  issues  involved. 

The  first  and  most  basic  simplification  has  usually 
already  been  applied  before  the  stability  problem  itself  has  even 
been  formulated.  It  consists  of  a  suitable  choice  of  the  slmpll- 
fied  energy  expression  (3»l«l)  on  which  the  ensuing  Investigation 
of  stability  is  based.  The  adequacy  of  the  chosen  energy  expres¬ 
sion,  or  of  an  equivalent  formulation  of  the  problem  in  terms  of 
equations  of  motion  and  stress-strain  re3ations,  is  often  taken 
for  granted.  Such  complacency  is  Indeed  more  often  than  not 
justified,  but  some  exceptions  prove  the  case  for  extreme  caution. 
It  should  for  Instance  always  be  borne  in  mind  that  actua.1 
problems  of  elastic  stability  concern  three-dimensional  elastic 
bodies,  no  matter  whether  one  or  two  dimensions  of  the  body  are 
small  in  comparison  with  the  other  dimension(s) .  The  representa¬ 
tion  of  the  potential  energy  of  the  body  by  a  sui’face  integral 
or  line  integral  in  the  case  of  approximately  two-dimensional  or 


92 


one -dimensional  structures  always  merits  a  comparison  with  the 
basic  three-dimensional  theory. 

A  sound  physical  insight  is  of  course  by  far  the  best 
guide  in  devising  appropriate  simplifications  of  the  basic  energy 
expression  (3^1.l).  V/ithout  such  insight  or  intuition  it  would 
be  virtually  hopeless  to  attempt  the  solution  of  any  problem  in 
elastic  stability.  Nevertheless^  it  should  be  kept  in  mind  that 
physical  Insight  is  no  more  than  a  precipitation  of  past  experl - 
ence,  and  it  may  well  become  unreliable;,  if  extrapolation  from 
this  experience  is  carried  too  far.  A  careful  mathematical 
evaluation  of  the  Implications  of  all  simplifying  physical 
assimiptlons  is  therefore  always  desirable ^  if  not  actually 
indispensable.  Complete  rigour  in  a  mathematical  sense  cannot 
usually  be  achieved  here^  but  every  effort  should  be  made  to 
establish  the  limitations  on  the  validity  of  the  theory  as  a 
consequence  of  the  assumptions  in  question. 

The  foregoing  more  or  less  philosophical  digressions 
may  again  be  illustrated  by  means  of  the  examples  which  have  been 
discussed  before.  The  energy  expression  (3.1.^)  is  the  exact 
equivalent  of  Euler’s  famous  formulation  and  investigation  of 
the  stability  of  thin  rods  in  axial  compression  [7].  We  shall 
discuss  this  problem  in  detail  in  chapters  and  .  It 

will  appear  that  Euler’s  theory  is  completely  justified,  provided 
that  the  cross-section  of  the  bar  is  compact .  On  the  other  hand; 
an  entirely  different  mode  of  instability  may  become  more 
critical  for  bars  v;ith  a  thin-v/alled  open  cross-section. 

The  phenomenon  of  buckling  of  such  bars  in  a  combined  flexural 
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and  torsional  mode  had  been  entirely  overlooked  until  less  than 
30  years  ago . 

A  second  instructive  lesson  to  be  gained  from  the  first 
example  concerns  the  structure  of  the  energy  functional.  In 
(3.1.4)  it  consists  of  one  term  which  represents  the  elastic 
strain  energy  in  bending ^  and  a  second  term  which  represents  the 
energy  of  the  external  loads .  The  elastic  strain  energy  due  to 
compression  of  the  bar  has  been  neglected.  The  second  variation 
of  the  energy  functional  (3.2.11)  consists  again  of  two  terms ^ 
the  first  one  representing  the  positive  definite  second  variation 
of  the  elastic  strain  energy^  the  second^  negative  term  represent¬ 
ing  the  second  variation  of  the  energy  due  to  the  external  loads . 
Several  misconceptions  in  more  general  buckling  problems^  in 
particular  for  flat  plates^  stem  from  unwarranted  generalizations 
of  the  experience  v/ith  the  buckling  of  bars.  First  of  all,  it 
was  conjectured  that  the  second  variation  of  the  elastic  energy 
is  always  positive  definite  and  that  the  possibility  of  instability 
would  therefore  arise  only  in  problems  v/here  the  second  variation 
of  the  energy  of  the  external  loads  is  negative  [4].  This  mis¬ 
conception  has  even  now  not  quite  disappeared.  Remnants  of  it 
are  still  to  be  found  in  current  textbooks  on  elastic  stability, 
in  spite  of  the  fundamental  researches  of  Trefftz  [18,19]  and 
Marquerre  [11,12]  in  the  thirties  of  the  present  century.  This 
misunderstanding  is  closely  related  to  a  second  common  erroneous 
conjecture  in  the  buckling  of  flat  plates,  where  it  is  often 
assumed,  explicitly  or  tacitly,  that  extensions  and  shears  in 


r- 
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the  middle  surface  of  the  plate  may  be  neglected  in  buckling^  as 
long  as  the  deflections  remain  (infinitely)  small.  The  connec¬ 
tion  between  these  misconceptions  becomes  clear;,  if  the  origin 
of  the  second  variation  of  the  energy  of  the  external  loads  in 
(3.2.11)  is  examined  more  closely.  It  stems  from  the  second 
term  in  (3.1.4) ^  which  is  itself  a  consequence  of  the  assumption 
that  the  center  line  is  inextensible .  If  this  assumption  is 
dropped j  the  second  variation  of  the  energy  of  the  loads  dis¬ 
appears  entirely^  and  it  is  replaced  by  an  additional  term  of 
equal  magnitude  in  the  second  variation  of  the  elastic  strain 
energy j  due  to  extension  of  the  center  line  (cf.  section  ). 
This  alternative  formulation  is  in  agreement^  as  it  should  be, 
with  cur  second  example  for  the  general  case  of  dead  loading  in 
which  the  second  variation  of  the  energy  (3.2.17)  is  entirely 
due  to  the  elastic  strain  energy. 

As  a  last  example  of  a  simplified  energy  expression 
(3.1.1)  vie  discuss  the  background  of  (3.1.8)  for  the  elastic  body 
under  dead  loading.  This  expression  represents  an  approximation 
of  the  energy  increase  for  an  elastic  body  which  obeys  a  partic¬ 
ular  generalization  of  Hookers  law  for  finite  deformations.  The 


approximation  is  here  embodied  in  the  quadratic  terms  in  the 
strain  components  Before  the  simplification  in  question, 

they  appear  in  the  form 


— V— I 
l-2v 


(3.9.1) 


IJhk 


where  all  components  of  the  tensor  AE, 


are  small  in  modulus 
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compared  with  the  sheax"*  modulus  G.  The  small  additional  terms  in 
(3.9 •!)  2-re  equivalent  to  small  changes  in  the  elastic  constants 
of  the  material,  which  are  not  known  with  great  precision  anyway. 
On  physical  gi'-ounds  we  have  therefore  little  hesitation  in  omit¬ 
ting  the  small  additional  terms .  A  more  mathematical  justifica¬ 
tion  is  also  readily  available  in  this  case.  Because  the  first 
term  in  (3.9.1)  is  positive  definite  it  is  indeed  ensured  that 
the  additional  terms  never  exceed  in  absolute  value  a  cei'^tain 
small  fraction  of  the  principal  term.  This  argument  would  fail, 
however,  if  the  first  tem  were  not  definite.  In  such  a  case 
the  neglection  of  the  additional  termis  would  lack  an  adequate 
mathematical  substantiation.  Even  though  the  coefficients 
AE^jhk  small  compared  with  the  coefficients  in  the  first 

term,  in  the  case  of  an  indefinite  or  semi-definite  first  term, 
the  additional  terms  might  well  be  comparable  with  or  even 
dominant  over  the  first  term  for  some  choice  of  the  strain 
components.  It  is  indeed  always  sound  practice  in  the  discussion 
of  "small"  terms  in  our  energy  expressions  to  compare  them  with 
a  positive  definite  aggregate  of  similar  terms  in  the  same 


variables ,  whenever  such  a  comparison  is  at  all  possible. 

Once  the  basic  simplified  energy  expression  (3.I.I) 
has  been  firmly  established,  the  next  step  in  the  investigation 
of  stability  deals  with  the  examination  of  the  second  variation 
?^[u],  A  further  simplification  of  this  second  variation  is 
again  desirable  in  most  cases,  and  it  may  Indeed  often  be  obtained 
without  detrimental  effects  on  the  accui'acy  of  the  predictions 
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on  stability.  The  need  of  proper  care^  hov;ever^  can  hardly  be 
overemphasized.  Physical  intuition  is  less  effective  here  than 
in  the  establishment  of  a  suitable  simplified  complete  energy 
expression  (3.I.I). 

As  an  example j  we  may  again  discuss  the  second  varia¬ 
tion  (3.2.17)  of  the  energy  in  the  case  of  the  three-dimensional 
elastic  body  under  dead  loading.  All  tenns  in  the  integrand 

depend  on  the  displacement  gradients  u.  .  as  arguments.  The 

1  >  J 

coefficients  in  the  first  term  are  certainly  all  small 

compared  v;ith  the  shear  modulus  G  appearing  as  a  factor  of  the 

_2 

second  term;  the  order  of  magnitude  of  is  at  most  10  to 

10  for  most  engineering  materials  in  the  elastic  range.  The 
entire  stability  problem  would  be  obliterated,  however,  if  the 
first  term  in  (3.2.17)  were  unsuspectingly  omitted.  The  remain¬ 
ing  expression  would  be  positive  definite,  and  no  possibility  of 
instability  would  arise  at  all.  It  is  evidently  impermissible 
here  to  neglect  the  supposedly  small  first  term.  The  reason  is 
of  course  that  the  second  term  in  the  integrand  is  not  a  positive 
definite  function  of  the  displacement  gradients ,  Omission  of 
the  first  term  can  therefore  not  be  supported  by  a  mathematical 
argument . 


This  becomes  even  more  clear,  if  we  introduce  in  addi¬ 
tion  to  the  (fictitious)  linear  strain  tensor,  defined  in  section 

(3.5) 


(3.9.2) 
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o 


the  (fictitious)  linear  rotation  tensor 


a 


ij 


(3.9.3) 


The  latter  skev;- symmetric  tensor  is  again  fictitious  as  a  rotation 
tensor  in  the  case  of  finite  displacement  gradients;  it  coincides 
with  the  rotation  tensor  only  in  the  case  of  an  infinitesimal 
displacement  field  u.  All  displacement  gradients  may  now  be 
expressed  in  the  linear  strain  and  rotation  tensors 


u 


i.  j 


(3.9.^^) 


and  the  second  variation  (3.2.17)  may  be  written  in  the  form 

P^Eu]  =  I  [|  s,j(6^^-Hc^^)(ej^+aj^)  +  ^  (ehh)^hdv. 

(3.9.5) 


It  is  now  perfectly  clear  that  the  first  term  in  (3.9.5) 
cannot  be  omitted  entirely  in  comparison  with  the  second  term. 

It  also  depends^  in  factj  on  the  rotation  tensor  (3.9*3)  which 
does  not  appear  in  the  second  term.  But  it  is  also  evident  from 

(3.9.5)  that  the  part 


0  0 
ij^ih^Jh 


(3.9.6) 


of  the  first  term  is  certainly  negligible  in  comparison  with  the 

positive  definite  second  term  in  the  linear  strain  components. 

Moreover^,  in  view  of  the  smallness  in  absolute  value  of  S.  ./G^ 

^  J 

it  may  also  be  argued  that  the  average  values  of  cc.  taken  in 

X  J 

a  mean  square  sense j  must  be  large  in  comparison  with  the  average 
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values  of  0^^  in  a  similar  sense^  if  the  possibility  of  instabil¬ 
ity  is  to  arise  at  all.  The  order  of  magnitude  of  this  ratio 
must  be  (G/S)^^  where  3  is  for  instance  the  largest  principal 
stress  in  absolute  value.  A  simplified  version  of  (3-9 *5)  Is 
then 

^2^-^  "  , 

(3.9.7) 


and  the  relative  error  Involved  appears  to  be  of  the  order  of 
1 

magnitude  (S/G)^.  The  actual  error  of  (3.9.6)  in  buckling 
problems  of  slender  structures  may  be  even  less  because  the 
neglected  terms 


ij 


+  64U.CI .ct.,  0 
ih  jh^  Ij  ih  jh 


(3.9.8) 


are  usually  of  the  same  order  of  magnitude  as  the  terms  (3.9.6). 
This  statement  will  be  substantiated  in  section  6.  It  is  based 
on  the  physical  fact  that  in  a  buckling  mode  the  linear  strain 
components  0^.^  which  appear  in  (3*9.8)  with  a  non-vanishing 
factor  are  small  compared  with  the  linear  strain  components 

which  dominate  in  the  second  term  of  (3.9*7)  and  determine  its 
order  of  magnitude.  The  actual  error  of  (3.9.7)>  in  comparison 
with  the  original  expression  (3.2.17)^  is  therefore  mostly  of 
order  of  magnitude  S/G  times  the  second  term.' 


in 

'The  foregoing  discussion  vrhlch  resulted  in  the  modified  form 
(3.9.7)  of  the  second  variation  should  not  be  taken  to  imply 
that  the  latter  form  would  always  be  the  more  convenient  one 
for  practical  applications.  The  original  expression  (3.2.17) 
is  often  more  appropriate^  for  example  in  the  discussion  of 
neutral  equilibrium  in  section  3*5. 
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We  shall  now  assume  that  we  have  obtained  an  expression 
for  the  second  variation  of  the  energy  in  the  simplest  possible 
form  for  the  investigation  of  the  stability  problem  at  hand.  Even 
then  it  will  usually  not  be  possible  to  discuss  this  second 
variation  by  means  of  a  rigorous  mathematical  analysis .  In  most 
cases  we  shall  have  to  resort  to  approximate  methods,  either  in 
a  direct  discussion  of  the  second  variation  or  in  the  solution 
of  the  associated  linear  and  homogeneous  differential  equations . 

We  shall  not  enter  here  into  a  comparison  of  the  numerous  avail¬ 
able  methods  .  Apart  from  a  few  words  of  caution  with  respect  to 
further  simplification  of  the  differential  equations,  we  shall 
be  content  with  a  brief  discussion  of  the  well-knov/n  Rayleigh- 
Ritz  method  for  an  approximate  direct  solution. 

Tlie  eigenvalue  problem  in  terms  of  the  differential 
equations  and  boundary  conditions  equivalent  to  the  variational 
equation  (3.^.1l)j  and  the  similar  equations  of  neutral  equilib¬ 
rium  equivalent  to  (3.5.3) ;»  often  Invite  to  a  further  simplifica¬ 
tion  by  the  neglection  of  certain  "small"  terms.  Particular 
vigilance  is  needed  before  it  is  even  considered  to  yield  to 
such  a  temptation.  It  is  admittedly  sometimes  justified  to 
neglect  certain  "small"  terms  in  these  equations.  In  most  cases, 
however,  this  is  due  to  the  fact  that  the  second  variation  of  the 
energy  has  not  yet  been  reduced  to  its  simplest  appropriate  form 
"Small"  terms  in  the  differential  equations  whose  omission  can¬ 
not  be  justified  by  a  further  permissible  reduction  of  the  second 
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variation  of  the  energy  should  alv;ays  be  regarded  with  great 
suspicion.  It  should  always  be  kept  in  mind  that  a  differential 
equation  requires  the  vanishing  of  a  sum  of  terms ^  and  a 
supposedly  small  additional  term  may  have  an  appreciable  effect 
on  the  solution.  The  discussion  of  small  terms  in  the  second 
variation  of  the  energy  is  on  much  safer  ground,  if  it  can  be 
based  on  a  comparison  with  a  positive  definite  group  of  similar 
terms.  Nevertheless,  it  can  often  not  be  avoided,  in  the  Interest 
of  further  progress  toward  an  approximate  solution,  to  neglect 
some  "small"  terms  in  the  differential  equations  without  adequate 
prior  justification.  In  such  cases  it  should  always  be  attempted 
to  obtain  a  Justification  a  posteriori.  This  may  sometimes  be 
done  simply  by  a  careful  examination  of  the  extent  to  which  the 
approximate  solution  actually  violates  the  original  equations  . 

In  other  cases  the  required  justification  may  be  obtained  from 
the  theory  of  asymptotic  integration  of  differential  equations 
containing  a  small  parameter. 

By  far  the  most  Important  tool  for  the  approximate 
solution  of  problems  in  elastic  stability  is  provided  by  the 
direct  Rayleigh-Rltz  method  of  the  calculus  of  variations.  This 
method  was  originally  devised  by  Rayleigh  [l4]  for  the  approxi¬ 
mate  solution  of  the  problem  of  the  fundamental  mode  of  free 
vibrations.  It  v;as  already  applied  with  conspicuous  success  by 
Timoshenko  to  the  solution  of  a  wide  variety  of  stability 
problems,  '  long  before  the  energy  theory  of  elastic  stability 
was  put  on  a  firm  foundation  by  Trefftz  [18,19]. 

’  The  reader  is  referred  to  Timoshenko's  treatise  on  elastic 
stability  [1?]  for  detailed  references  to  his  numerous  papers. 
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The  Rayleigh -Rltz  method  for  the  Investigation  of  the 
second  variation  is  based  on  the  minimuin  problem  (3. '4. 2)  or 
(3.^.3).  In  the  latter  form  it  is,  for  an  appropriate  choice  of 
TgLuJj  even  identical  to  Rayleigh  *s  original  f ormiilation .  The 
minimum  property  of  the  (normalized)  minimizing  displacement  field 
is  characterized  by  (cf.  section  3.4) 


+  O(s^) 


(3.9.9) 


where  ^  is  an  arbitrary  kinematically  admissible  displacement 

o . 

field,  and  0(e‘"j  stands  for  a  positive  term  which  tends  to  zero 
2 

as  e  ,  if  6  approaches  zero.  The  Ray 1 eight -Ritz  method  now 
consists  of  the  evaluation  of 


w 


_  ^2^^^ 
“  THuT 


(3.9.10) 


for  some  suitably  selected  kinematically  admissible  displacement 
field,  which  is  Intended  to  represent  a  reasonable  approximation 
to  the  minimizing  displacement  field  .  The  result  of  the 

- X 

evaluation  of  (3.9.10)  is  taken  to  provide  an  approximation  of  the 
minimum  . 

Tlie  extraordinary  power  of  the  Rayleigh-Ritz  method 
lies  in  the  minimum  property  expressed  by  (3.9.9).  It  implies 
that  a  first-order  error  in  the  guess  for  the  minimizing  displace- 

t 

ment  field  results  in  a  second-order  error  in  .  Moreover,  more 
severe  local  errors  in  the  guessed  displacement  field  are  evened 
out  by  the  integrations  involved  in  both  numerator  and  denominator 

of  (3.9.10). 
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It  follows  iminediatsly  from  (3.9*9)  that  the  approximate 
value  always  exceeds  the  actual  minimimi  .  The  approximate 
verdict  on  stability  obtained  by  the  RaylelpihrRitz  method  always 
overestimates  the  stability.  If  equilibri^mI  is  unstable  accord- 

I 

ing  to  the  approximate  verdict that  is  if  Wt  is  negative^  this 
equilibrium  is  certainly  actually  mstable.  On  the  other  hand, 

I 

a  positive  value  of  is  no  guarantee  for  actual  stability. 

It  needs  hardly  to  be  emphasized  that  a  skillful  choice 
for  the  displacement  field  u  in  (3.9.IO)  is  of  paramount  import¬ 
ance  for  the  achievement  of  an  accurate  prediction  on  stability. 
Considerable  scope  exists  here  again  for  the  deployment  of  past 
experience  and  physical  insight.  Some  assistance  may  also  be 
obtained  from  the  differential  equations  and  dynamic  boundary 
conditions j  even  if  they  have  defeated  attempts  at  a  complete 
rigorous  solution.  In  particular  it  is  often  not  too  difficult 
to  choose  the  displacement  field  u  in  (3.9.10)  in  such  a  way 
that  it  satisfies  the  dynamic  boundary  conditions,  in  addition 
to  the  geometric  boundary  conditions .  The  latter  conditions 
must  of  course  always  be  satisfied  by  a  kinematically  admissible 
displacement  field. 

The  chances  of  a  satisfactory  result  are  furthermore 
greatly  enhanced,  if  we  incorporate  in  the  displacement  field  u 
a  number  of  suitably  selected  indeterminate  parameters  c^^ 
(k=l,2,...).  Upon  substitution  of  the  displacement  fields 
into  (3.9.10),  the  best  possible  approximation  is  obtained  as 
the  solution  of  the  algebraic  minimum  problem 
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“1  ■ 

The  expression  to  be  minimized  is  here  a  function  of  the 
parameters  The  labour  involved  increases  very  rapidly  with 

the  number  of  parameters,  and  a  good  physical  judgment  is  again 
Invaluable,  if  satisfactory  results  are  to  be  obtained  without 
excessive  effort.  In  Ritz*s  special  form  of  the  method  the 
parameters  c,  all  enter  into  the  displacement  field  u(c,  )  in  a 
linear  manner.  This  field  is  then  described  by 

u  =  2  c.  ,  (3.9.12) 

k  ^ 

(k)  / 

where  each  of  the  chosen  displacement  fielcis  u^  '  (which  should 
not  be  confused  v/lth  the  eigenfunctions)  satisfies  the  geometric 
boundary  conditions .  This  form  of  the  method  has  the  advantage 
that  both  numerator  and  denominator  in  (3.9*ll)  are  homogeneous 
quadratic  functions  in  the  parameters  c^^.  The  resulting  condi¬ 
tions  for  a  minimum  of  ( 3.9.11)  are  then  homogeneous  and  linear 
algebraic  equations,  for  which  many  effective  methods  of  solution 
are  available.  Nevertheless,  it  is  often  more  convenient  to 
include  one  or  more  parameters  in  ^(cj^)  in  a  nonlinear  form,  for 
example  as  an  indeterminate  wave  length,  in  order  to  achieve  an 
accurate  approximation  of  the  minimizing  displacement  field  u^ 
with  the  smallest  possible  number  of  indeterminate  parameters. 

The  Rayleigh-Rit z  method  has  a  serious  drawback  in  that 
it  furnishes  only  an  upper  bound  for  the  minimum  .  In  order  that 
the  accuracy  of  an  approximation  may  be  established  with  certainty. 
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a  lower  bound  for  is  also  required.  Such  a  lower  bound  is 
even  more  Important  than  an  upper  bound  from  the  point  of  view 
of  safety  in  engineering  applications  because  it  permits  a 
conservative  prediction  on  the  stability  of  a  structure . 
Unfortunately,  it  is  a  far  more  difficult  problem  to  establish 
a  lower  bound  for  a  minimum  with  some  accuracy  [8,l6].  Vie 
shall  not  discuss  an^^  of  the  methods  devised  for  this  purpose 
because  their  application  is  hardly  worthv/hile  in  most  problems. 
Any  available  additional  effort  for  the  establishment  of  a  lower 
bound  is  usually  more  rewardingly  spent  on  an  improvement  of 
previously  established  upper  bounds . 

The  application  of  the  Rayleigh-Ritz  method  to  the 
determination  of  critical  loads  in  the  theory  of  elastic  stability 
is  conventionally  formulated  in  a  slightly  different  manner.  A 
critical  load  is  defined  here  as  an  external  load  system  under 
which  the  fundamental  state  is  a  critical  case  of  neutral 
equilibrium.  The  semi-definite  character  of  the  second  varia¬ 
tion  P^Lu],  resulting  in  a  solution  w^=0  of  minimum  problem 
(3.4.2),  may  alternatively  be  expressed  by  the  statement  that 
the  minimum  value  zero  of  P2[u]  is  attained  for  a  non-vanishing 
displacement  field  u^ .  In  the  case  of  our  first  example  of  the 
bar  under  end  loads,  whose  second  variation  of  the  energy  is 
given  by  (3.2.11),  the  critical  load  may  therefore  be 
characterized  by 
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N 


1 


Min. 


^  Bw 

1  ‘P 
^  w 


(3.9.13) 


with  geometric  boundary  conditions  w(0)=w(i)=0. 

A  similar  formulation  is  often  given  for  the  general 
three-dimensional  elastic  body  under  dead  loading^  whose  second 
variation  is  given  by  (3.2.1?).  The  initial  stresses  are 
then  specified  In  terms  of  a  unit  system  of  initial  stresses 
and  a  load  factor  \ 

Si^.  =xs[5) 


The  condition  for  a  semi-definite  second  variation  of  the  energy 
(3.2.1?)  may  then  be  written  in  the  form 

Mln.{J  + 

+  \  I  °  ’  (3.9.15) 

where  X=X^  denotes  the  critical  load  factor.  The  minimum  problem 
(3.9.1^)  is  nov;-  conventionally  replaced  by 


=  Min. 


-  /  i  sj^^u,  .u,  .dv 

^2  ij  h,l  h, j 


.  (3.9.16) 


Some  comments  are  called  for  in  connection  with  the 
traditional  formulation  of  the  critical  load  problem  in  the 
form  (3.9.16).  It  rests  on  a  first  tacit  assumption  that 
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(3.9 is  an  appropriate  representation  for  the  initial 
stresses.  This  assumption  is  Justified,  if  and  only  if  the 
change  in  geometry  of  the  body  in  the  fundamental  state  I, 
compared  with  the  undeformed  state,  may  be  neglected  (cf. 
chapter  6) .  In  such  a  case  we  have  indeed  (approximate)  pro¬ 
portionality  of  the  initial  stresses  with  the  load  factor,  and 
the  differences  in  geometry  of  the  body  in  state  I  for  different 
values  of  the  load  factor  may  be  neglected  in  the  evaluation  of 
the  integrals  in  (3.9.15)  and  (3.9.16).  The  case  in  v/hich  the 
fundamental  state  is  a  state  of  pure  strain  (cf.  chapter  6) 
represents  an  important  group  of  examples  in  v;hich  this  assump¬ 
tion  is  Justified.  In  fact,  it  is  Justified  in  most  of  the  more 
important  problems  of  elastic  stability,  but  exceptions  do  exist, 
and  it  should  always  be  verified  whether  the  application  of 
(3.9.15)  or  (3.9.16)  is  legitimate  in  any  particular  case. 

A  second  tacit  assumption  is  made  in  the  transition 
from  (3.9.15)  to  (3.9.16).  It  is  assumed  here  that  the  denomina¬ 
tor  in  (3.9.16)  is  positive  definite,  and  this  assimiption  is 
Incorrect  in  a  wide  class  of  problems .  It  often  happens  that 
both  positive  and  negative  critical  load  factors  e^xist,  in  which 
cases  the  denominator  in  (3. 9-16)  is  indefinite.  No  solution 
exists  of  problem  (3.9.16)  in  such  cases  because  the  functional 
in  (3.9.16)  is  unbounded.  The  obvious  remedy  here  is  to  consider 
Instead  of  (3.9.16)  the  reciprocal  maximum  and  minimum  problems 
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1/X^  =  Max. 
1A“  =  Min. 


(3.9.17) 


If  the  body  Is  properly  supported,  the  denominator  in  (3.9.1?) 
is  positive  for  all  non-vanishing  displacement  fields .  The 
maximLim  value  of  (3.9.17)  Is  then  the  reciprocal  value  of  the 
positive  critical  load  factor  and  the  minimum  of  (3-9.17)  is 
the  reciprocal  value  of  the  negative  critical  load  factor  X“. 
Equilibrium  in  the  fundamental  state  is  stable  for  X“  <  X  <  X^  . 

The  Raylelgh-Ritz  method  may  also  be  applied  to  the 


Itivestigation  of  a  critical  case  of  neutral  equilibrium. 


quantity  (3.6.40),  which  gives  the  verdict  on  stability  in 
the  case  of  a  simple  buckling  mode  u^,  depends  on  the  solution 
Vg  of  the  minimum  problem 

Pglv]  +  P2i[ui>v]  =  Min.  (3.9.18) 


under  the  side  condition 


=  0  .  (3.9.19) 

If  we  assume  a  displacement  field  v  which  satisfies  (3.9.19)> 
and  v/hich  contains  a  number  of  indeterminate  parameters,  we 
obtain  the  best  possible  approximation  v^^  if  we  minimize  (3.9.18) 
with  respect  to  these  parameters .  In  viev/  of  the  minimum 
property  of  the  exact  solution,  we  have  in 
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A4  =  P4U1]  -  Pgtvg]  (3.9.20) 

an  upper  bound  for  Here  £igain  the  Rayleigh-Ritz  method 

overestimates  stability. 

The  foregoing  discussion  presupposes  that  the  critical 
case  of  neutral  equilibrium  has  been  established  rigorously^ 
i.e.  that  the  exact  solution  oj^of  minimum  problem  (3.^.2)  is 
zero,  and  that  that  minimizing  displacement  field  u^  is  completely 
known.  If  we  have  already  previously  applied  the  approximate 
Rayleigh-Ritz  method  to  establish  an  approximation  for  the 
critical  load,  we  may  again  apply  this  method  in  order  to 
investigate  the  stability  of  the  (approximate)  critical  case  of 
neutral  equilibrium.  The  buckling  mode  employed  in  this  investi¬ 
gation  is  then  also  approximate.  Hence  our  investigation  of  the 
minimum  problem  (3*9«18)j  (3.9 *19) >  where  the  exact  buckling 
mode  has  been  replaced  by  its  approximation,  is  also  app.roxlmate 
in  the  sense  that  even  its  exact  solution  would  provide  no  more 
than  an  estimate  for  A|^.  Moreover,  we  then  no  longer  knov/ 
whether  this  estimate  is  too  large  or  too  small. 
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